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shifters  were  measured  and  statistically  analyzed.  Several  research  contracts  for  promising  phased  array  techniques 
were  awarded  to  industrial  contractors.  In  addition  to  research  on  the  primary  MLS  functions,  which  involve  scan¬ 
ning  linear  phased  array  antennas,  NRL  has  studied  the  realization  of  a  360-degree  scanning  antenna  utilizing  a  cir¬ 
cular  array.  This  report  contains  a  nontechnical  review  of  all  activities,  supported  by  technical  appendices  and  a 
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NRL  PARTICIPATION  IN  THE 
MICROWAVE  LANDING  SYSTEM  PROGRAM 


I .  INTRODUCTION 

1.1  P  UXg.P.S.6. 

This  report  summarizes  work  performed  by  the  Radar  Division, 
Naval  Research  Laboratory,  under  Interagency  Agreement  No.  DOT-FA75WAI- 
556.  This  agreement  was  initiated  in  May  1975  between  the  Federal 
Aviation  Administration  (FAA)  and  the  Naval  Research  Laboratory  (NRL). 
This  report  covers  activity  during  the  period.  May  1975  through  Septem¬ 
ber  1979. 

This  report  is  the  final  summary  report  (FY-79)  called  for  in  para¬ 
graph  L,  Article  V,  in  Modification  10  to  the  Interagency  Agreement. 

The  Interagency  Agreement  has  been  further  modified  to  extend  the  effort 
through  FY-80.  Therefore,  this  research  program  is  continuing,  and  this 
report  is  not  final  in  the  strict  sense.  However,  several  phases  of  the 
program  have  been  completed,  and  it  is  feasible  to  present  a  cohesive 
summary  of  activity  and  progress  to  date. 

The  general  objective  of  this  research  program  is  to  study  scanning 
beam  antenna  designs  and  to  make  recommendations  to  the  FAA  as  to  the 
applicability  of  these  designs  to  the  Microwave  Landing  System  (MLS). 

NRL  is  participating  in  this  program  because  of  its  expertise  in  the 
technology  being  utilized  in  MLS  antennas.  During  the  period  covered  by 
this  report,  the  MLS  program  devoted  very  little  of  its  resources  to 
military  applications  of  MLS,  and  NRL  has  devoted  only  a  very  slight 
fraction  of  its  effort  to  consideration  of  military  application  of  MLS. 
Thus,  this  research  program  has  been  concerned  with  MLS  as  an  advanced 
precision  landing  system  for  civil  aviation.  Appendix  A  contains  the 
Statement  of  Work  from  the  original  Interagency  Agreement  and  subsequent 
modifications  thereto. 

Although  the  MLS  program  is  predominantly  directed  to  civil  aviation 
applications,  DoD  and  the  Navy  clearly  have  strong  interest  both  in  the 
MLS  program  in  particular  and  in  advanced  air  traffic  control  technology 
in  general.  Thus,  it  is  considered  to  be  in  the  mutual  interest  of  FAA 
and  the  Navy  to  have  NRL  participate  in  this  program.  From  the  stand¬ 
point  of  the  Radar  Division,  where  the  NRL  work  on  this  program  has  been 
performed,  the  technology  interchange  relating  to  system  concepts  and 
detailed  design  techniques  has  been  valuable. 
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For  those  readers  who  are  not  familiar  with  the  Microwave  Landing 
System,  Appendix  B  is  included  as  a  brief  introduction  with  both  his¬ 
torical  and  technical  information. 

Unfortunately,  there  has  been  relatively  little  treatment  of  the 
MLS  in  the  open  literature.  The  background  and  general  guidelines  for 
the  U.S.  development  effort  are  given  in  the  final  report  of  RTCA  Spe¬ 
cial  Committee  117  [1],  The  FAA  sponsored  MLS  symposia  in  1972  and 
1973,  for  which  records  are  published  [2,3].  Some  general  descriptions 
of  landing  systems  can  be  found  [4,5,6].  Workers  from  Bendix  and 
Hazeltine  have  described  the  results  of  some  of  their  efforts  [7,8,9]. 
Articles  in  news-type  publications  such  as  Aviation  Week  and  Electronic 
News ,  or  even  the  general  news  media,  emphasized  the  adversary  aspects 
of  the  ICAO  process.  Other  articles  contained  self-admitted  bias  [10] 
or  were  based  on  a  general  misunderstanding  of  the  issues  and  processes 
involved  [11].  The  general  subject  of  air  traffic  control  is  receiving 
increased  popular  coverage,  and  the  MLS  program  was  recently  mentioned 
in  National  Geographic  Magazine  [12]. 

1.2  Outline  of  NRL  Activity 

Because  of  NRL's  participation  in  the  MLS  program  in  the 
capacity  of  expert  consultants,  because  of  NRL's  close  proximity  to  FAA, 
and  because  of  the  wide  range  of  activities  being  pursued  by  FAA  during 
the  time  period  covered  by  this  report,  a  variety  of  tasks  were  under¬ 
taken  by  NRL  in  accomplishing  the  objectives  of  the  statement  of  work. 

As  has  been  previously  stated,  the  initial  research  objective  was  to 
make  recommendations  to  FAA  as  to  the  applicability  of  scanning  antenna 
designs  to  the  MLS  program. 

At  the  outset,  NRL  addressed  microwave  optical  techniques  for  the 
reasons  that  (a)  the  Phase  III  programs  being  undertaken  by  Bendix  and 
Texas  Instruments  (TI)  were  both  using  microwave  optical  designs  and 
(b)  the  Australian  INTERSCAN  MLS,  with  which  the  U.S.  had  aligned  it¬ 
self  at  the  conclusion  of  the  selection  process  at  the  end  of  1974, 
used  microwave  optical  designs.  Therefore,  intensive  interaction  was 
required  between  NRL  and  Bendix,  TI,  and  the  Radio  Physics  Laboratory 
of  the  Commonwealth  Scientific  and  Industrial  Research  Organization  of 
Australia  (CSIRO).  Part  of  this  interaction  was  achieved  through  the 
medium  of  the  Design  Review  Board,  to  which  NRL  was  a  consultant,  which 
reviewed  the  progress  of  Bendix  and  TI  on  the  Phase  III  program  during 
1975-76. 

A  basic  purpose  for  making  antenna  recommendations  to  FAA  was  to 
strengthen  the  U.S.'s  submission  to  ICAO,  which  was  initially  made  in 
December  1975.  Further,  the  ICAO  submission  was  defended  and  modified 
during  the  ICAO  meetings  of  1976-77. 

NRL  evaluated  proprietary  proposals  from  ITT  Gilfillan  and  Hazel- 
tine  in  1975  and  eventually  awarded  research  contracts  to  those  companies 
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for  studies  of  novel  phased  array  techniques.  These  studies  were  com¬ 
pleted  in  1976,  and  a  second  contract  was  awarded  to  Hazeltine  for  mod¬ 
ification  and  additional  testing  of  its  technique.  A  study  contract  was 
awarded  to  Myer  Associates  in  1979  to  assist  with  NRL's  performance  and 
monitoring  simulations. 

As  a  result  of  various  cost  and  technical  factors,  emphasis  shifted 
from  microwave  optical  techniques  to  phased  arrays  in  early  1976.  Con¬ 
sequently,  NRL  also  devoted  increasing  effort  to  phased  array  studies, 
initially  considering  the  design  and  performance  of  thinned  arrays,  and 
later  addressing  the  monitoring  and  performance  simulation  of  all  types 
of  digitally  scanned  arrays. 

1.3  Introduction  to  Report 

This  report  is  organized  with  the  hope  of  providing  the  reader 
with  a  broad  overview  of  NRL's  participation  in  the  MLS  program  during 
the  1975-77  time  period.  Although  no  attempt  is  made  to  completely  sat¬ 
isfy  all  types  of  readers,  it  is  hoped  that  the  information  contained 
herein  will  be  useful  to  a  wide  variety  of  readers,  ranging  from  engi¬ 
neers  with  specific  design  interests  to  those  individuals  with  more 
general  interest  in  the  recent  history  and  evolution  of  the  MLS  program. 

Sections  3  and  4  are  a  glossary  of  terms  and  references,  respec¬ 
tively.  Sections  5  and  6  list  documents  generated  by  NRL  and  personnel 
participating  in  the  program.  Appendices  A,  B,  and  C  include  the  work 
statement  from  the  interagency  agreement,  a  brief  description  of  the 
Microwave  Landing  System,  and  a  description  of  contract  research  sup¬ 
ported  by  NRL. 

Appendix  D  is  a  listing  in  chronological  order  of  the  significant 
events  and  milestones  relating  to  NRL's  activities  in  the  MLS  program. 

In  many  cases  the  milestones  are  also  highly  significant  to  the  overall 
program,  so  that  the  chronology  reasonably  reflects  the  U.S.  and  ICAO 
MLS  activity  during  the  period. 

Appendix  E  is  a  description  of  the  performance  simulation  which  has 
been  developed  by  Dr.  J.  K.  Hsiao,  and  Appendix  F  is  an  analysis  of  the 
digital  filter  used  in  that  simulation. 

Appendix  G  describes  the  optimum  beam  steering  technique  generated 
at  NRL.  Appendix  H  described  the  circular  array  design  study  carried  out 
in  1978-79.  Appendix  I  describes  3D  bootlace  lenses  which  were  developed 
on  another  NRL  task  but  which  have  direct  application  to  the  original 
microwave  optical  scanning  studies  and  the  circular  array  study.  Appen¬ 
dix  J  describes  how  to  feed  a  linear  array  with  a  3D  lens. 

2.  DISCUSSION  OF  ACTIVITIES 

2,1  Microwave  Optical  Research 
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The  original  statement  of  work,  agreed  upon  between  NRL  and  FAA  and 
reproduced  in  Appendix  A,  listed  the  following  seven  subjects  for  con¬ 
sideration. 


1.  Microwave  optics. 

2.  Feed  commutation. 

3.  Monitoring  and  redundancy. 

4.  Beam  shaping. 

5.  Phased-array /microwave  optics  trade-off  studies. 

6.  Alternate  proposal  for  ICAO. 

7.  Design  review  activities. 

The  emphasis  was  very  strongly  on  scanning  antennas  utilizing  microwave 
optics.  The  Phase  III  contracts  were  just  underway  in  mid-1975,  and 
both  TI  and  Bendix  were  using  microwave  optical  designs  for  the  Basic 
Narrow  (narrow  referring  to  the  aperture  size)  configuration,  in  which 
the  azimuth  antenna  beamwidth  is  about  two  degrees  and  the  elevation 
antenna  beamwidth  is  about  1.5  degrees.  Bendix  used  Gent  bootlace 
lenses  (also  known  as  Rotman)  for  their  Phase  III  equipment,  and  TI 
used  a  bifocal  folded  pillbox  for  the  elevation  unit  and  a  parabolic 
torus  for  the  azimuth  unit. 

Although  these  microwave  optical  designs  were  felt  to  be  cost 
effective  in  the  Basic  Narrow  configuration,  they  were  not  considered 
practical  for  the  Basic  Wide  configuration,  which  requires  3.66m  (12- 
foot)  apertures  to  prod.uce  one-degree  beamwidths.  Thus,  two  objectives 
of  the  microwave  optical  studies  were  (a)  to  review  all  possible  micro- 
wave  optical  techniques  to  determine  what  would  be  the  preferred  design 
approach,  and  (b)  to  ascertain  whether  there  was  a  practical  microwave 
optical  design  for  use  in  the  Basic  Wide  configuration. 

In  addition  to  the  microwave  optical  design  activity  underway  at 
Bendix  and  TI  in  1975,  the  Australian  INTERSCAN  MLS  program  was  committed 
to  microwave  optical  scanning  antennas.  The  principal  MLS  research 
activity  in  Australia  has  taken  place  at  the  Radio  Physics  Laboratory 
of  CSIRO  in  Sydney,  under  the  leadership  of  Dr.  J.  P.  Wild  and,  later, 

H.  Minnett.  The  principal  MLS  design  and  fabrication  contractor  has 
been  Australasian  Wireless  Amalgamated. 

Immediately  upon  execution  of  the  interagency  agreement,  it  was 
deemed  desirable  to  coordinate  the  research  programs  at  NRL  and  CSIRO. 

To  this  end,  meetings  were  held  during  a  two-week  period  in  June  1975 
in  Paris  and  Sydney.  The  American  representatives  in  these  discussions 
were  J.  P.  Shelton  of  NRL  and  S.  R.  Jones  of  MITRE.*  Initial  and  final 
discussions  were  held  with  J.  P.  Wild  in  Paris,  and  detailed  technical 

discussions  were  held  with  H.  Minnett,  associate  director  of  the  Radio 

Physics  Laboratory,  and  other  workers.  The  outcome  of  these  meetings 

was  an  agreement  to  pursue  an  intensive  cooperative  antenna  study  pro¬ 

gram  during  the  period  July-September  1975,  with  another  meeting  at  the 

*M [TRK  corporation  was  an  overall  technical  consultant  to  FAA  on  the 
MLS  program,  and  S.  R.  Jones  was  the  MITRE  project  leader  until  1977. 


4 


end  of  the  study  to  decide  upon  an  appropriate  hardware  development 
program.  The  Australians  described  two  antenna  concepts  which  they 
regarded  as  proprietary:  One  was  the  use  of  curved  waveguide  slot 
array  radiators  in  circular  azimuth  arrays,  which  eliminates  beamwidth 
broadening  at  high  elevation  angles.  The  other  was  a  cylindrical  geo¬ 
desic  dielectric-loaded  lens,  which  appeared  equivalent  to  the  Myers 
geodesic.  The  agreed-upon  study  program  included  the  following  objec¬ 
tives  : 

1.  Antennas  with  one-degree  beamwidth,  obtained  with  an  effective 
aperture  of  75A.  In  the  case  of  circular  configurations,  the  feed 
radius  will  be  no  greater  than  40A. 

2.  CSIRO  would  pursue  (a)  complex  modulation  analysis  and  simula¬ 
tion,  (b)  experimental  improvement  of  torus  antenna  to  demonstrate  one- 
degree  beamwidth,  (c)  design  of  column  waveguide  arrays  and  shaped 
reflectors,  and  (d)  study  of  new  disk  and  cylinder  lenses. 

3.  NRL  would  pursue  (a)  general  lens  studies,  (b)  feed  network 
concepts,  and  would  provide  any  possible  information  on  lens  references 
and  slot  array  design. 

4.  Both  groups  agreed  to  study  monitoring,  redundancy,  cost,  and 
error  budget. 

It  was  estimated  that  the  next  meeting  would  be  about  mid-September. 
The  subsequent  meeting  actually  occurred  on  1  December  1975,  at  CSIRO 
in  Sydney.  The  American  group  consisted  of  J.  P.  Shelton,  NHL,  S.  R. 
Jones,  MITRE,  and  G.  Jensen  and  R.  Normoyle  of  FAA.  This  two-week 
meeting  produced  significant  results,  although  they  were  not  predicted 
by  the  objectives  listed  above.  The  cooperative  antenna  study  program 
was  remarkable  in  producing  innovative  concepts;  it  is  doubtful  that 
these  innovations  would  have  been  generated  without  the  impetus  of  the 
cooperative  program. 

The  following  developments  were  presented  by  CSIRO: 

1.  J.  P.  Wild  described  his  "Opera  House"  lens,  so  named  because 
of  its  resemblance  to  the  Sydney  opera  house.  This  bootlace  lens  con¬ 
figuration  utilizes  a  spherical  geodesic  surface  rather  than  the  conven¬ 
tional  planar  internal  propagation  region.  It  has  been  shown  that 
Wild's  lens  is  the  optimum  trifocal  bootlace  configuration. 

2.  Torus  efforts  included  experimental  development  of  an  eight- 
phase  complex  modulator,  plus  theoretical  design  of  a  compact  torus 
with  R  *  77. 5\  (the  objective  was  75A). 

3.  Parallel-plate  lens  designs  for  feeding  circular  arrays  in¬ 
cluded  flat  disk  lenses  and  spherical  geodesic  configurations  with 
linear  angle  transformations. 


The  following  developments  were  presented  by  NRL: 
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1.  Redundant  switching  and  modulation  networks  had  been  developed 
in  mid-October  for  inclusion  in  the  ICAO  submission.  Both  types  of  net¬ 
work  rely  on  PIN  diode  control  devices  configured  so  that  the  network 
can  sustain  the  failure  of  one  diode  without  hindering  operation.  The 
switches  use  hybrid  matrix  networks  with  typically  4,  8,  or  16  ports  on 
each  side.  Available  switch  configurations  for  these  matrices  are  SP3T, 

SP7T,  and  SP15T  with  nonisolated  outputs  or  SP2T,  SP4T,  and  SP8T  with 
all  ports  isolated  and  matched.*  The  redundant  modulator  technique  uses 
diode  phase  shifters  with  a  Fourier  transforming  Butler  matrix  or  equiv¬ 
alent  lens  to  transform  from  phase  scan  to  lateral  feed  motion. 

2.  A  simplified  analysis  of  Gent-Rotman  lenses  was  carried  out  to  j 

establish  the  performance  limitations  of  the  lens.  It  was  concluded  that  ] 

the  Gent-Rotman  lens  can  be  used  to  generate  one-degree  beams  with  +40  j 

degrees  coverage  with  the  required  precision.  It  was  also  concluded 

that  the  optimum  lens  design  has  front-back  symmetry. 

3.  A  modular  Rotman  lens  configuration  was  described,  which  uses 
small  lenses  to  build  larger  lenses  analogous  to  the  building-block 
assembly  of  Butler  matrices. 

4.  Parallel-plate  lenses  with  circular  symmetry  were  described. 

Technical  memoranda  on  (a)  a  family  of  folded  disk  lenses,  (b)  a  geo¬ 
desic  Luneburg  design,  and  (c)  a  bootlace-corrected  trifocal  Schmidt 
system  were  presented.  To  aid  in  the  realization  of  parallel-plate  lens 
structures,  a  parallel-plate  directional  coupler  was  proposed. 

At  the  conclusion  of  the  December  1975  meeting  at  CSIRO,  plans  were 
made  to  progress  from  Stage  I  of  the  cooperative  antenna  program  into 
Stage  II,  which  would  involve  experimental  testing  of  some  of  the  pro¬ 
posed  novel  techniques. 

Although  work  on  microwave  optical  techniques  continued  at  CSIRO, 
effort  at  NRL  was  diverted  from  microwave  optics  into  phased-array  tech¬ 
nology  early  in  1976.  Subsequent  efforts  are  described  in  the  following 
sections. 

2.2  Phased-Array  Studies 

2.2.1  Background  for  Phased-Array  Effort 

The  cooperative  NRL-CSIRO  microwave  optical  study  effort 
produced  notable  results.  However,  in  terms  of  the  broader  requirements 


*It  was  subsequently  learned  that  this  switch  configuration  was  not  orig¬ 
inal.  The  basic  matrix  switching  configuration  was  proposed  by  H.  Schrank 
in  the  60's  (Schrank,  Hooper,  and  David,  "A  Study  of  Array  Beam  Switching 
Techniques,"  RADC-TR-65-546 ,  Feb.  65)  and  C.  Lee  more  recently  reinvented 
the  SP3T  version  (U.S.  Pat.  3,996,533,  7  Dec  76). 
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of  the  U.S.  MLS  program  and  the  presssures  of  the  ICAO  selection  process, 
it  was  decided  in  early  1976  to  adopt  phased  arrays  as  the  preferred 
design  approach  for  the  U.S.  TRSB  MLS.  The  following  factors  contribu¬ 
ted  to  this  decision: 

1.  Although  techniques  for  achieving  redundancy  in  the  switching 
and  modulation  portions  of  the  microwave  optical  antenna  systems  were 
proposed,  it  became  clear  that  experimental  validation  and  field  testing 
of  these  techniques  could  not  be  accomplished  soon  enough  to  make  an 
impact  on  the  ICAO  process, 

2.  The  cost  of  phase  shifters  was  dropping  dramatically,  with  the 
result  that  the  cost  advantage  of  microwave  optical  systems  was  becoming 
much  less  significant. 

3.  Although  microwave  optical  antennas  could  be  designed  to  meet  the 
MLS  requirements  for  one-degree  beamwidths,  the  lens  structures  were 
always  larger  than  the  aperture  length.  For  the  one-degree  beamwidth, 
this  leads  to  integral  planar  structures  with  dimensions  of  about  3.66m 
(12  feet),  which  is  excessively  bulky  for  equipment  which  must  be  trans¬ 
portable  on  a  world-wide  basis. 

NRL’s  involvement  in  these  phased-arrav  studies  can  oe  divided 
roughly  into  three  stages.  The  first  stage  consisted  of  a  theoretical 
investigation  of  thinned  arrays.  This  work  coincided  with  experimental 
field  tests  of  the  Eendix  azimuth  phased  array  at  N’AFEC,  modified  to  a 
chinned  configuration.  The  second  stage  of  the  URL  phased-array  effort 
consisted  of  a  study  of  integral  RF  monitoring  techniques  for  fault 
detection  and  isolation,  carried  out  through  computer  simulation.  The 
third  stage  involved  the  development  of  a  performance  simulation  appli¬ 
cable  to  all  of  the  phased  array  configurations  under  consideration.  In 
addition,  NRL  has  established  a  procedure  for  testing  phase  shifters 
that  are  intended  for  use  in  MLS  phased  arrays. 

2.2.2  Thinned  Arrays 


It  was  decided  early  in  1976  to  prepare  a  case  for  the 
use  of  thinned  phased  arrays  at  airports  where  multipath  is  not  a  serious 
problem.  The  thinned  array  is  realized  by  quasi-random  removal  of  ele¬ 
ments  from  a  filled  array.  The  resultant  system  can  be  made  to  exhibit 
essentially  unchanged  pattern  shape  near  the  main  beam  at  the  cost  of 
much  higher  far-out  sidelobes,  together  with  the  lower  gain  associated 
with  a  smaller  number  of  radiators. 


The  questions  to  be  answered  during  the  thinned  array  study  related 
to  the  effects  of  thinning  on  accuracy,  reliability,  multipath  vulnera¬ 
bility,  monitorability ,  and  cost.  These  topics  were  treated  in  a  report 
titled  "Preferred  Time  Reference  Scanning  Beam  Phased  Array  Implementa¬ 
tion  for  MLS,"byS.R.  Jones,  MITRE,  J.P.  Shelton,  NRL,  and  R.E.  Willey, 


Bendix.  The  report  was  entered  as  a  Background  Information  Paper  to  the 
Washington  meeting  of  Working  Group  A  of  AWOP,  held  in  May  1976.  NRL 
contributed  about  two-thirds  of  this  report  with  the  following  sections: 

2.0  Design  Concepts 

2.1  Aperture  Features  and  Beam  Scanning  Effects 

2.1.1  Sidelobe  Performance  of  Thinned  Arrays 

2.1.2  Beam  Scanning  Effects  in  Phased  Arrays 

2. 1.2.1  Effect  of  Phase  Shifter  Bit  Reduction  on  Sidelobe  Level 

2. 1.2.2  The  Beam-Pointing  Accuracy  of  a  Phased  Array 

2.2  Failure  Characteristics  of  Phased  Arrays 

3.0  Hardware  Considerations 

3.2  Advanced  Techniques 

3.2.1  Integrated  Diode-Controlled  Phased  Array  Technique 

3.2.2  Network  Techniques  for  Limited  Scan  Arrays  (COMPACT) 

Appendix  A:  A  Review  of  Thinning  Techniques  for  Linear  Arrays 
(by  J.  K.  Hsiao) 

Meanwhile,  preparations  were  underway  at  Bendix  to  perform  accuracy 
measurements  of  a  thinned  azimuth  array  at  NAFEC.  The  array  was  the  one 
developed  by  Bendix  during  Phase  II  in  the  1973/74  time  period.  Thinning 
was  accomplished  by  removing  phase  shifters  from  the  filled  array.  The 
performance  of  the  thinned  array  had  been  predicted  by  simulations  carried 
out  at  Bendix.  The  measured  performance,  obtained  in  a  series  of  tests 
carried  out  in  June  1976,  was  in  good  agreement  with  the  simulation 
results.  Furthermore,  phase  shifter  failures  were  measured  and  simulated, 
and  the  graceful  degradation  of  the  phased  array,  even  in  the  thinned 
condition,  was  amply  demonstrated.  These  test  results  were  submitted  to 
the  AWOP  WGA  meeting  held  at  The  Hague  in  July  1976. 

2.2.3  Phased  Array  Monitoring  Studies 

After  the  July  WGA  meeting  at  The  Hague,  the  U.S.  TRSB  pro¬ 
posal  was  modified  to  include  both  thinned  and  filled  phased  arrays,  and 
emphasis  was  placed  on  analysis  of  a  fault-locating  monitoring  system  for 
phased  arrays. 

A  monitoring  simulation  program  was  begun  at  NRL  in  August  1976. 

The  obiective  was  to  complete  the  simulation  so  that  the  results  could  be 
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used  at  the  London  WGA  meeting  scheduled  for  November  1976.  The  simula¬ 
tion  was  sufficiently  complete  so  that  the  fault  isolation  capability 
of  the  proposed  U.S.  monitoring  system  could  be  confirmed.  The  outcome 
of  the  London  meeting  was  the  desired  one:  that  monitoring  was  not  a 
discriminant  between  DHLS  and  TRSB  MLS. 

2.2.4  Phased  Array  Performance  Simulation 

After  the  AWOP  meeting  at  Montreal  in  May  1977  and  the  vote 
by  AWOP  in  favor  of  the  U.S.  system,  a  decision  was  made  by  the  MLS 
office  to  develop  an  "in-house"  performance  simulation  to  provide  the 
capability  for  comparing  the  various  phased  array  configurations  and 
beam-scannin2  algorithms  being  proposed  by  industry.  It  was  decided  that 
NRL  would  develop  this  simulation,  and  effort  was  begun  in  April  1977. 
Appendix  E  describes  the  present  status  of  this  simulation. 

2.2.5  Phase  Shifter  Measurements 

In  the  fall  of  1977  FAA  became  interested  in  the  performance 
that  could  be  expected  from  phase  shifters  that  might  be  used  in  MLS 
antennas  and  the  effect  of  component  tolerances  on  system  performance. 

In  1975  Bendix  had  purchased  180  phase  shifters  from  EMS  for  possible 
use  in  the  Phase  III  equipment  development  program.  Although  these  units 
were  not  used  in  Phase  III  equipment,  it  was  decided  that  they  would  be 
used  in  the  Basic  Wide  system. 

Therefore,  NRL  was  tasked  to  measure  a  small  number  of  phase  shift¬ 
ers,  about  10,  ac  the  outset,  then  a  group  of  about  30,  and  finally  the 
entire  available  inventory  of  180. 

Although  these  units  were  purchased  under  a  set  of  specifications, 
it  was  necessary  to  evaluate  their  actual  performance  so  that  the  mea¬ 
sured  performance  of  the  system  could  be  related  to  the  actual  phase 
shifter  characteristics. 

Initially,  a  group  of  about  seven  Bendix/ EMS  ferrite  latching  four- 
bit  units  and  two  TI  diode  four-bit  units  were  tested  using  an  automatic 
network  analyzer  available  at  NRL.  A  computer  program  was  written  for 
compiling  statistical  data  on  the  measurements,  which  reads  out  inser 
sertion  phase  characteristics  and  standard  deviation  of  the  variation  of 
phase  shift  from  the  desired  values,  together  with  mean  and  standard 
deviation  of  insertion  loss  and  VSWR.  Subsequently,  NRL  was  tasked  with 
testing  130  Bendix/EMS  units.  At  this  point  it  was  realized  that  the 
manual  transfer  of  data,  by  means  of  punched  cards,  from  the  printed 
output  of  the  network  analyzer  to  the  computer  was  laborious,  time- 
consuming,  and  error-prone.  Therefore,  a  somewhat  circuitous  but  fully 
mechanized  procedure  was  devised,  whereby  the  test  data  was  recorded  on 
magnetic  tape,  then  transcribed  to  a  disk  and  input  to  the  computer. 
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Testing  and  analysis  of  the  180  units  was  performed  in  FY-79.  It  was 
planned  that  a  subset  of  this  group  would  be  retained  at  NRL  for  tem¬ 
perature  testing.  Unfortunately,  Bendix  called  for  return  of  these 
remaining  units  before  the  testing  could  be  completed. 

2.2.6  Optimum  Beam  Steering  Technique 

In  the  course  of  the  research  effort  on  phased  arrays,  it  was 
noted  that  all  contractors  had  proprietary  beam  steering  algorithms. 
Subsequently,  J.  P.  Shelton  of  NRL  proposed  an  optimum  beam  steering 
technique  which  uniquely  determines  the  sequence  in  which  phase  shifters 
will  be  stepped.  This  technique  is  described  in  Appendix  G. 

2 . 3  Investigation  of  360-degree  Ground  Scanning  Azimuth  Antenna 

Following  the  selection  by  ICAO  in  Ap^il  1973  of  the  TRSB/MLS 
and  a  decision  at  that  meeting  that  the  design  of  a  360-degree  ground 
scanning  azimuth  antenna  should  be  developed,  a  portion  of  NRL's  effort 
has  been  devoted  to  this  problem. 

In  'lay  1978  J.  P.  Shelton  of  NRL  proposed  a  three-dimensional  boot¬ 
lace  lens  configuration  which  is  applicable  to  the  generation  of  multiple 
beams  from  linear  arrays.  This  development  has  impact  on  two  problem 
areas  in  the  MLS.  It  had  already  been  found  that  some  of  the  feed  sys¬ 
tems  which  would  be  used  for  the  360-degree  system  require  lenses  with 
dimensions  larger  than  those  of  the  circular  array.  The  3D  lens  will 
have  dimensions  much  smaller  than  those  of  the  corresponding  2D  lens. 

Thus,  it  is  possible  to  use  a  3D  lens  in  the  360-degree  antenna  system 
which  is  smaller  than  the  circular  array. 

Furthermore,  the  3D  lens  solves  the  original  problem  which  NRL  was 
tasked  with  in  1975.  It  is  now  possible  to  feed  the  one-degree-beam- 
width  azimuth  array  with  a  lens  structure  which  is  small  relative  to  the 
size  of  the  array. 

Appendix  H  is  a  discussion  of  the  circular  array  design  study. 
Appendix  I  describes  two  3D  lens  configurations  which  have  been  derived. 
Appendix  J  gives  the  procedure  for  applying  the  3D  lens  to  a  linear 
array. 

2.4  Design  Review  Board  and  Contractor  Interaction 

In  order  to  monitor  the  progress  of  the  Phase  III  development 
contracts  at  Bendix  and  TI,  the  MLS  office  organized  a  Design  Review 
Board  (DRB)  which  was  active  during  the  critical  phases  of  the  contracts — 
from  Che  summer  of  1975  through  the  spring  of  1976.  Meetings  were  held 
at  the  contractors'  various  facilities  at  roughly  monthly  intervals. 

NRL  participated  in  the  DRB  process  as  an  expert  consultant  and  provided 
an  attendee  to  most  of  the  DRB  meetings. 
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In  addition,  NRL  acted  as  a  source  of  consulting  expertise  to  the 
contractors.  Indeed,  the  interaction  between  NRL  and  the  contractors 
was  of  great  benefit  to  all  parties.  The  degree  of  interaction  was  more 
or  less,  depending  on  such  factors  as  location  or  the  disposition  of  the 
contractor  to  discuss  his  activities.  (Some  contractors  show  a  marked 
disinclination  to  discuss  technical  matters  except  when  a  contract  is  in 
the  offing.)  The  results  of  these  factors  were  as  follows  for  the  var¬ 
ious  contractors: 

Bendix  -  Strong  because  of  both  location  and  disposition  and  also 
heavy  and  continuing  MLS  activity. 

TI  -  Weak  primarily  because  of  location. 

ITT  Gilfillan  -  Weak  because  of  remote  location  and  brief  period 
of  activity. 

Hazeltine  -  Negligible  because  of  disposition. 

2.5  ICAO  Participation 

It  has  already  been  pointed  out  that  NRL  submitted  material  for 
inclusion  in  the  U.S.  MLS  ICAO  proposal,  which  was  submitted  in  December 
1975.  To  some  extent,  all  MLS  activity  by  FAA  has  been  directed  to  the 
ICAO  process  because  of  the  absolute  necessity  of  winning  that  contest. 
Subsequent  to  the  submission  of  proposals  to  ICAO  by  Australia,  UK,  FRG, 
and  US,  there  were  four  AWOP  WGA  meetings  in  1976  followed  by  an  AWOP 
meeting  in  March  1977. 

NRL  participated  in  the  preparation  of  the  thinned  phased  array  pro¬ 
posals  that  were  entered  into  the  ICAO  process  at  the  Washington  meeting 
in  May  1976  and  at  The  Hague  meeting  in  July  1976. 

J.  P.  Shel~on  of  NRL  attended  the  London  WGA  meeting  in  November 
1976  and  the  Montreal  AWOP  meeting  in  March  1977. 


Russian  delegations  visited  Washington  in  November  1975  and  November 
1977,  the  Russians  supported  TRSB  MLS.  NRL  participated  in  the  discussions 
with  the  Russian  delegations. 

2.6  Contracts  with  ITT  Gilfillan  and  Hazeltine 

In  mid-1975  when  NRL  entered  the  MLS  program,  the  Phase  III 
contracts  with  TI  and  Bendix  were  also  underway.  Thus,  all  of  the 
developmental  effort  being  funded  by  FAA  at  that  time  was  directed 
toward  microwave  optical  scanning  antennas,  and  there  was  no  effort 
being  directed  to  phased  arrays. 

At  roughly  the  same  time,  ITT  Gilfillan  and  Hazeltine  came  for¬ 
ward  with  unsolicited  proprietary  proposals  for  novel  phased-array 
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techniques.  It  was  decided  that  it  would  be  advantageous  to  the  MLS 
program  to  investigate  these  techniques.  NRL  was  deemed  the  logical 
organization  to  issue  and  monitor  the  necessary  contracts. 


One  technique,  proposed  by  ITT  Gilfillan,  employs  a  "vari-phase 
coupler"  to  integrate  the  power  distribution  to  the  array  elements  with 
a  diode  phase-shift  procedure.  The  technique  offers  the  advantages  of  a 
compact  physical  package,  low  losses,  and  low  cost  with  the  possible 
disadvantages  of  three-bit  phase  control  and  amplitude  variation  with 
phase  state.  This  integrated  phased  array  technique  is  described  in 
Appendix  C. 

ITT  Gilfillan  fabricated  a  44-element  horizontal  array  with  three- 
degree  beamwidth  under  their  contract.  The  array  is  fed  at  the  center. 

The  other  technique,  under  development  at  Hazeltine  and  acronymed 
COMPACT,  employs  a  network  between  the  array  and  the  phase  shifters  to 
form  near-optimum  sector  element  patterns  so  that  a  minimum  number  of 
phase  shifters  can  be  used  in  limited  scan  applications,  such  as  the 
MLS  elevation  antennas.  The  network  requires  that  the  number  of  phase 
shifters  be  one-half  the  number  of  inputs  to  the  array,  and  the  tech¬ 
nique  is  limited  to  coverage  of  about  40  to  50  degrees.  The  limited- 
scan  network  technique  is  also  described  in  Appendix  C. 

Hazeltine  fabricated,  under  their  research  contract,  a  96-element 
vertical  array  with  one-degree  beamwidth  which  can  be  scanned  by  24 
phase  shifters.  The  96  elements  consist  of  48  element  pairs,  which  are 
fed  by  the  COMPACT  network. 

NRL  issued  and  monitored  contracts  for  both  of  these  techniques, 
and  both  contractors  completed  their  contracts  in  mid-1976.  Final 
reports  were  submitted,  and  the  experimental  measurements  showed  excel¬ 
lent  agreement  with  the  performance  predicted  for  the  techniques.  In 
addition,  NRL  awarded  a  second  contract  to  Hazeltine  in  April  1977  for 
additional  testing  of  the  COMPACT  antenna. 

In  July  1979  NRL  issued  a  contract  to  Meyer  Associates  for  a  simu¬ 
lation  and  monitoring  study.  This  study  encompassed  the  following  tasks 

1.  Review  the  NRL  MLS  performance  simulation  program  and  develop 
procedures  for  validating  it. 

2.  Study  the  performance  of  the  field  monitor  under  conditions  of 
component  failures  and  develop  a  suitable  simulation  procedure  and  vali¬ 
dation  process. 

3.  Review  present  and  future  requirements  for  maintenance  monitor¬ 
ing  in  the  MLS  ground  scanning  antennas. 

4.  Design  an  experimental  maintenance  monitor  for  evaluating  phase 
shifter  performance. 
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3.  GLOSSARY  OF  TERMS 


ALPA 

ATA 

AWA 

AWOP 

COMPACT 

CSIRO 

DLS 

DMLS 

DRB 

FAA 

FRSB 

IATA 

ICAO 

IFALPA 

ILS 

INTERSCAN 

MLS 

NAFEC 


NRL 

TDM 

TI 

TRSB 

WGA 


Airline  Pilots'  Association 

Air  Transport  Association 

Australiasian  Wireless  Amalgamated 

All-weather  Operations  Panel 

Cost  Minimized  Phased  Array  Circuit  Technique 
(proprietary  to  Hazeltine) 

Commonwealth  Scientific  and  Industrial  Research 
Organization  (Australia) 

DME-based  Landing  System  (West  Germany) 

Doppler  Microwave  Landing  System  (U.K.) 

Design  Review  Board 

Federal  Aviation  Administration 

Frequency  Reference  Scanning  Beam 

International  Air  Transport  Association 

International  Civil  Aviation  Organization 

International  Federation  of  Airline  Pilots'  Associations 

Instrument  Landing  System 

Australian  TRSB  MLS 

Microwave  Landing  System 

National  Airport  Facilities  Evaluation  Center, 

Atlantic  City,  N.J. 

Naval  Research  Laboratory 

Time  Division  Multiplex 

Texas  Instruments 

Time  Reference  Scanning  Beam 

Working  Group  A 
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APPENDIX  A 


Interagency  Agreement  Work  Statement 
and  Modifications  Thereto 

The  original  work  statement  contained  in  the  Interagency  Agree¬ 
ment  of  May  1975  is  as  follows: 

III.  STATEMENT  OF  WORK 

NRL  shall  provide  all  the  required  personnel,  facilities  equipment, 
materials  and  services  necessary  to  perform  the  effort  herewith  set 
forth: 

A.  This  development  effort  shall  include  assessments  of  antenna 
designs  and  performance,  consideration  of  alternative  design  techniques, 
and  trade-off  analyses. 

1.  Microwave  Optics.  Focusing  systems  for  microwave  optical  scan¬ 
ners  shall  be  investigated.  Both  wide-  and  narrow-angle  scan  sectors 
shall  be  considered,  as  well  as  linear  and  circular  apertures.  An 
extensive  investigation  shall  be  made  of  scanning  microwave  optical  sys¬ 
tems,  and  their  applicability  to  the  MLS  shall  be  assessed.  Where  nec¬ 
essary,  as  in  the  case  of  the  torus,  for  example,  a  ray  tracking  computer 
program  shall  be  generated  to  assess  pattern  performance.  An  attempt 
shall  be  made  to  establish  the  general  limitations  of  circularly  sym¬ 
metric  collimators,  with  a  view  toward  relating  the  number  of  degrees  of 
freedom  available  in  the  structure  to  the  wavefront  tolerance,  and  thereby 
the  limiting  aperture  size, 

2.  Feed  Commutation.  The  Technology  of  beam  scan  by  commutating  a' 
feed  array  shall  be  reviewed.  Relationships  among  beam  shape  continuity 
and  linearity  of  scan,  feed  radiator  density,  feed  array  length  and 
mutual  coupling  effects  shall  be  examined.  The  possibility  of  aberra¬ 
tion  correction  by  commutation  power  division  shall  be  investigated. 
Methods  of  minimizing  insertion  losses  shall  be  sought, 

3.  Monitoring  and  Redundancy.  Internal  monitoring  systems  shall 
be  studied  for  all  scanning  techniques  of  interest.  Monitoring  systems 
shall  be  assessed  in  terms  of  accuracy,  complexity,  and  reliability. 

The  details  of  location  of  monitoring  components  shall  be  worked  out  for 
the  different  types  of  antenna  configurations.  Redundancy  shall  be 
studied  from  the  standpoints  of  failure  modes  of  commutated  feed  systems 
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and  required  parallel  stand-by  components.  The  effect  of  redundant 
design  on  cost  and  performance  shall  be  analyzed. 

Beam  Shaping.  The  effect  of  the  requirement  for  elevation 
beam  shaping  on  the  design  of  the  azimuth  antenna  shall  be  analyzed, 

A  ray-tracing  computer  program  shall  be  used  to  compute  patterns  for 
circular  arrays  of  column  radiators,  for  a  linear  aperture  feeding  a 
cylindrical  shaped-beam  reflector,  and  for  a  toroidal  aperture  (see 
III. A. 1. ) . 

5.  Phased-Array/Microwave  Optics  Trade-off  Studies.  This  trade¬ 
off  analysis  shall  evaluate  a  number  of  commutating  techniques  and 
phased-array  designs  according  to  several  criteria,  such  as  size,  ini¬ 
tial  cost,  operating  cost,  logistics,  reliability,  monitorability ,  and 
maintainability . 

6.  Alternate  Antenna  Proposal  for  ICAO.  NRL  shall  prepare  an 
alternate  antenna  proposal  to  be  included  in  the  U.S.  MLS  submission 
to  ICAO.  This  proposal  shall  specify  prime  alternate  antenna  designs 
appropriate  for  TRSB  MLS;  a  preferred  alternate  design  approach  shall 
be  identified  taking  into  account  the  results  of  other  antenna  design 
activity  in  the  MLS  Program.  This  choice  shall  be  well  supported  by 
available  analyses,  simulations,  bench  tests  and  in-depth  considera¬ 
tions  of  monitoring,  redundancy  requirements,  cost  factors  and  possible 
failure  modes.  In  sum,  this  effort  shall  yield  a  well  thought  out  pro¬ 
posal.  to  ICAO  of  desirable  alternate  TRSB  antenna  designs. 

7.  Design  Review  Activities.  NRL  shall  participate  in  Design 
Review  activities  as  appropriate.  These  reviews  will  be  held  with  the 
two  prime  industrial  teams  at  3-6  week  intervals  throughout  the  MLS 
phase  III.  The  FAA  anticipates  NRL  attendance  at  those  Design  Reviews 
where  antenna  configurations  are  a  primary  topic.  This  activity  shall 
include  comment  and  critique  of  these  configurations. 

In  September  1975  the  following  paragraph  was  added  to  the  Work 
Statement  by  means  of  Modification  No.  1: 

8.  Analysis  of  New  Antenna  Concepts.  NRL  shall  evaluate  new 
antenna  designs  showing  promise  for  cost  savings  in  phased  arrays. 

The  purpose  of  this  modification  was  to  provide  for  NRL's  activity  in 
phased  array  studies  and  in  supporting  outside  contracts  for  promising 
phased  array  techniques. 

In  February  1977  with  Modification  No.  5,  the  Work  Statement  was 
further  modified  as  follows: 


17 


lo  Paragraph  A. 2.,  after  the  last  sentence  add  the  following: 

"Study  beam  steering  algorithms  and  required  logic  hardware 
for  phased  arrays  including  fully  filled,  limited  scan,  such 
as  the  Cost  Minimized  Phased  Array  Circuit  Technique  (COMPACT), 
and  thinned  configurations." 

2.  Paragraph  A. 3.,  after  the  last  sentence  add  the  following: 

"Conduct  studies  of  rf  monitoring  techniques  for  phased  arrays. 
Address  system  requirements,  analyze  and  simulate  proposed 
techniques,  and  fabricate  breadboard  systems  as  required. 
Establish  appropriate  rf  component  tolerances  to  match  the 
operation  of  the  monitor  to  the  required  operating  performance 
of  the  Time  Reference  Scanning  Beam  (TRSB)  antenna." 

3.  Paragraph  A.80,  as  contained  in  Modification  1,  after  the  last 
sentence  add  the  following: 

"Antennas  acquired  shall  be  tested  to  compare  the  performance 
with  the  low  cost  design  goals.  Review  the  current  state-of- 
the-art  in  phase  shifter  technology." 

In  October  1977  with  Modification  No.  7  the  Work  Statement  was 
further  modified  as  follows: 

1.  Paragraph  A2  as  amended  by  Modification  No.  5,  after  the  last 
sentence  add  the  following: 

"Digital  simulation  of  phased  array  antenna  configurations 
will  be  utilized  to  investigate  current  MLS  designs.  NRL  will 
provide  comment  and  critique  to  FAA  with  respect  to  system 
design  considerations." 

2.  Add  the  following  Paragraph  A9. 

"NRL  shall  test  components  such  as  slotted  wave  guide  radiators, 
phase  shifters,  and  shall  perform  microwave  network  analysis. 
Data  shall  be  taken  for  inclusion  in  a  report." 

In  March  1979  with  Modification  No.  10  the  Work  Statement  was  fur¬ 
ther  modified  as  follows: 

Article  IV.  Deliverable  Items,  add  the  following  paragraphs 

K.  Draft  Summary  Report  -  (FY-79) 

Submit  three  copies  of  the  draft  report  documenting  the 
results  of  the  work  accomplished  during  the  period  October 
1,  1978  through  September  30,  1979,  on  Items  III  A  3  and  3. 
Prepare  in  NRL  format. 
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L.  Final  Summary  Report  -  (FY-79) 

Submit  one  reproducible  copy  of  the  report  generated  under 
the  Task  in  Item  III  A  2  and  3. 

Article  V.  Completion  of  Work  &  Delivery,  Paragraph  I,  Draft  of 
Network  Analysis  Report,  as  amended  in  Modification  No.  7 

"Delete  July  31,  1978,  and  in  lieu  therefor  insert  March 
31,  1979" 

Paragraph  J.  Network  Analysis  Report,  as  amended  in  Modification 
No.  7 

"Delete  September  30,  1978  and  in  lieu  thereof  insert 
May  31,  1979." 

Add  the  following  paragraphs: 

"K.  Draft  Summary  Report  (FY-79) 

Submit  not  later  than  October  30,  1979.  The  FAA  will 
require  thirty  (30)  days  to  review,  approve,  disapprove, 
or  request  changes  to  this  report." 

"L.  Final  Summary  Report  (FY-79) 

Submit  not  later  than  December  31,  1979." 
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APPENDIX  B 


A  Brief  Description  of  the  Microwave  Landing  System 

History  -  The  first  commercial  ILS  was  demonstrated  in  1939,  The 
use  of  ILS  expanded  rapidly,  and  it  was  adopted  by  ICAO  as  the  inter¬ 
national  standard  in  1949.  Efforts  to  develop  a  second-generation 
landing  guidance  system  were  initiated  in  the  late  1950's.  As  a  result 
of  a  request  by  ATA  to  the  FAA  in  1967,  special  committee  117  (SC-117) 
of  the  Radio  Technical  Commission  for  Aeronautics  was  formed  in  December 
1967  to  develop  a  precision  guidance  system  concept  for  approach  and 
landing.  The  U.S.  Government  formed  an  interagency  planning  group  to 
establish  a  coordinated  national  plan  for  the  development  of  the  Micro- 
wave  Landing  System.  The  resulting  plan  provided  for  FAA  leadership 
in  developing  an  MLS  suitable  for  submission  to  ICAO  as  an  international 
standard. 

In  1973  FAA  funded  the  development  of  doppler  and  scanning  beam 
systems.  In  December  1974  Time  Reference  Scanning  Beam  was  selected  as 
the  United  States  candidate  MLS.  TRSB  MLS  was  submitted  to  ICAO  by  the 
FAA  in  December  1975.  The  All  Weather  Operations  Panel  (AWOP)  of  ICAO 
voted  in  favor  of  TRSB  MLS  in  March  1977,  and  ICAO  ratified  the  AWOP 
decision  at  its  world  wide  meeting  in  April  1978. 

Thus,  NRL  has  been  participating  in  the  MLS  program  under  this 
Interagency  Agreement  since  shortly  after  the  TRSB  selection  was  made. 
Some  further  details  of  MLS  programs  both  in  the  U.S.  and  elsewhere 
during  the  time  period  from  about  1973  to  the  present  are  germane  to 
this  discussion. 

The  FAA  funding  of  the  development  of  both  doppler  and  scanning 
beam  systems  was  referred  to  as  Phase  II  of  the  MLS  program.  The  dop¬ 
pler  systems  were  designed  by  ITT  Gilfillan  and  Hazeltine,  The  scanning 
beam  systems  were  designed  by  Texas  Instruments  and  Bendix.  The  scan¬ 
ning  beam  systems  utilized  frequency  reference  (FRSB),  and  TI  used 
mechanically  scanning  antennas,  while  Bendix  used  electronically  scanned 
phased  arrays.  These  systems  were  installed  and  tested  at  NAFEC  and 
Wallops  Island. 

Meanwhile,  MLS  development  was  being  actively  pursued  by  other  gov¬ 
ernments.  Australia  was  developing  an  electronic  scan  TRSB  system 
referred  to  as  INTERSCAN.  The  U.K.  has  pioneered  doppler  MLS  (DMLS) 
and  was  developing  that  system.  France,  West  Germany,  and  Russia  were 
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developing  systems,  but  with  somewhat  less  intensity.  West  Germany 
evenutally  proposed  their  DMR-based  landing  system  (DLS)  to  ICAO. 


This  then  was  the  situation  in  the  fall  of  1974  when  the  FAA 
entered  its  selection  process,  which  included  representatives  from  other 
agencies  of  the  Government,  such  as  NASA  and  DoD,  from  other  governments 
such  as  Australia,  U.K.,  France,  and  West  Germany,  and  from  other 
organizations  such  as  ATA  and  ALPA. 

It  was  realized  during  the  selection  process  that  TRSB  is  superior 
to  FRSB  from  the  standpoint  of  spectrum  occupancy,  although  TRSB  pre¬ 
cludes  the  use  of  mechanically  scanned  antennas.  The  option  of  mechani¬ 
cal  scan  was  dropped,  and  FRSB  was  replaced  by  TRSB  as  the  scanning 
beam  candidate.  Finally,  TRSB  was  selected  over  doppler. 

Subsequent  to  the  selection  process  the  FAA  contracted  with  Bendix 
and  TI  in  1975  for  Phase  III  systems,  which  involved  electronically 
scanned  microwave  optical  antennas  and  which  were  derived  from  the 
Australian  INTERSCAN  system. 

Finally,  as  the  cost  of  microwave  electronic  phase  shifters  designed 
specifically  for  the  MLS  application  fell  dramatically,  the  FAA  placed 
increasing  emphasis  on  phased  arrays,  especially  for  the  most  capable 
systems,  beginning  in  1976. 

The  selection  of  TRSB  by  the  U.S.  was  viewed  with  alarm  and  con¬ 
sternation'  by  the  U.K.  There  had  apparently  been  a  basic  misunderstand¬ 
ing  relative  to  the  FAA  selection  process:  The  U.S.  assumed  that  the 
U.K.  would  align  its  position  to  the  candidate  selected  by  the  FAA,  and 
the  U.K.  assumed  that  the  U.S,  would  select  the  doppler  system. 

Neither  the  U.S.  nor  the  U.K.  apparently  considered  the  alternative 
outcomes,  and  the  result  of  this  shortsightedness  was  a  bitterly  fought 
contest  for  adoption  by  ICAO,  which  was  waged  through  the  four  Working 
Group  A  meetings  of  1976,  through  the  AWOP  meeting  of  March  1977  (at 
which  the  vote  was  in  favor  of  TRSB),  and  continued  unabated  into  the 
world  wide  ICAO  meeting  which  was  held  in  April  1978.  At  that  meeting 
the  U.S./Australian  TRSB/MLS  was  the  winner. 

Technical  -  A  technical  description  of  MLS  is  best  introduced  by 
stating  the  general  requirements  of  the  system.  The  highest  category 
of  service  to  be  provided  by  MLS  is  IIIC,  as  defined  by  ICAO.  This  will 
provide  guidance  information  to  the  aircraft  of  sufficient  precision  to 
allow  automatic  landing.  MLS  provides  azimuth  angle  information  over  a 
region  +60  deg.  relative  to  runway  center  line  with  vertical  coverage 
to  +20  deg.  Elevation  angle  information  is  provided  from  0  to  20  deg. 
Back  azimuth  or  missed  approach  information  is  provided  over  a  region 
+40  deg.  relative  to  runway  center  line.  Flare  elevation  information 
is  provided  from  -2  to  +8  deg.  Data  is  also  provided  to  identify  func¬ 
tion,  airport,  runway,  nominal  glide  slope,  minimum  glide  slope,  and  the 
like. 
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The  US,  and  the  Australian  TRSB  and  the  U.K.  DMLS  are  closely 
related  in  several  respects.  The  operating  frequency  is  C-band,  5030 
to  5090  MHz.  Transmitters  are  on  the  ground,  and  receivers  are  in  the 
aircraft;  angle  information  is  derived  in  the  airborne  receiver  and 
displayed  to  the  pilot. 

The  basic  difference  between  the  doppler  and  TRSB  systems  is  that 
the  doppler  MLS  codes  the  antenna  in  time  and  codes  the  coverage  region 
in  frequency,  while  TRSB  codes  the  antenna  in  frequency  and  codes  the 
coverage  region  in  time.  One  is  essentially  the  Fourier  transform  of 
the  other,  or  in  other  words  one  is  the  dual  of  the  other  in  terms  of 
time  and  frequency,  DMLS  measures  frequency  in  its  airborne  receiver, 
and  TRSB  measures  time.  For  the  same  antenna  apertures,  both  systems 
will  have  basically  the  same  capability. 

Therefore,  the  differences  between  the  doppler  and  scanning  beam 
techniques  are  found  in  the  second-order  design  characteristics,  such 
as  the  requirement  for  doppler  to  use  two  transmitters,  one  as  a  refer¬ 
ence,  or  differences  in  data  rates  and  signal  formats. 

To  elaborate  slightly  on  the  operation  of  the  scanning  beam  system, 
the  most  capable  version  will  use  four  scanning  antennas,  each  of  which 
generates  a  fan  beam.  The  azimuth  and  back  azimuth  antennas  will  produce 
vertical  fans  which  scan  to  and  fro  in  azimuth  to  provide  azimuth 
information  to  an  airborne  receiver.  The  elevation  and  flare  antennas 
will  have  horizontal  fan  beams  which  scan  to  and  fro  in  elevation  to 
provide  elevation  information  to  an  airborne  receiver.  Suggested  bean- 
widths  for  the  fan  beams  are  one  degree  for  azimuth  and  elevation,  three 
degrees  for  back  azimuth,  and  one-half  degree  for  flare.  The  resultant 

aperture  sizes  are  roughly  3.6  m  (12  feet),  1.2  m  (4  feet),  and  7,3  m 
(24  feet),  respectively.  The  azimuth  antenna  is  located  on  center  line 
at  the  far  end  of  the  runway,  and  the  back  azimuth  is  located  at  the 
side  of  the  runway  opposite  the  desired  touchdown  region.  The  flare 
antenna  is  located  at  the  side  of  the  runway  some  distance  beyond  the 
touchdown  region. 

Angle  information  is  time  division  multiplexed  (TDM),  and  the 
transmitted  rf  signal  is  switched  sequentially  among  the  various 
antennas.  For  each  scanning  antenna,  the  beam  is  scanned  in  one  direc¬ 
tion  and  then  the  other.  The  time  difference  between  the  arrivals  of 
the  beam  maxima  of  the  to  and  fro  scans  is  measured  in  the  airborne 
receiver  and  is  directly  related  to  the  angular  position  of  the  air¬ 
craft.  Each  scanning  beam  transmission  is  identified  by  a  preceding 
data  transmission. 

Low  cost  has  been  one  of  the  primary  objectives  of  MLS  antennas, 
because  the  system  will  be  in  widespread  commercial  use  and  it  is 
planned  that  one  version  will  be  applicable  to  small  community  airports. 
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High  precision  is  of  course  a  dominant  requirement.  Scanning 
antenna  accuracy  of  the  order  of  1/50  beamwidth  is  required  to  provide 
the  specified  overall  system  precision. 

Because  the  system  is  involved  with  the  most  critical  portion  of 
the  approach  and  landing  procedure,  it  must  not  introduce  any  hazard. 
Thus,  reliability,  monitorability ,  maintainability,  redundancy,  and 
failure  tolerance  are  all  subjects  of  crucial  importance  to  the  program. 

Finally,  the  concept  of  maintainability  is  important  in  two 
respects.  First,  the  system  must  be  easily  and  quickly  maintainable 
so  that  it  has  a  minimum  of  down  time.  Second,  it  must  be  easily  main¬ 
tainable  by  relatively  unskilled  personnel  because  it  will  be  deployed 
world  wide,  and  there  must  be  no  possibility  of  system  malfunction  be¬ 
cause  of  improper  maintenance. 


\ 
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APPENDIX  C 

Contractor  Research  Supported  by  NRL 

This  appendix  contains  qualitative  descriptions  of  the  concepts 
proposed  by  ITT  Gilfillan  and  Hazeltine  and  supported  by  NRL  contracts. 
Both  of  the  techniques  are  novel,  and  the  emphasis  in  this  appendix  is 
on  a  clear  presentation  of  the  basic  concept. 

1.  E-PLANE  VARIPHASE  EXCITER 

A  diagram  of  an  E-plane  variphase  exciter  is  shown  in  Figure  C-l. 

The  device  employs  deep  slots  in  the  narrow  wall  of  a  waveguide  that 
intercept  the  longitudinal  currents  in  the  main  series  guide.  When  the 
diode  pair  in  a  slot  is  symmetrically  biased  in  the  forward  or  reverse 
state,  the  net  coupling  to  the  radiating  element  is  zero.  This  is  so 
because  a  slot  intercepts  equal  and  opposite  rf  currents  on  the  top  and 
bottom  walls  of  the  main  guide,  which  results  in  a  net  coupling  of  zero 
to  the  TE^g  mode  of  the  output  guide.  If  now  one  diode  of  the  pair  is 
reverse  biased,  and  the  other  diode  is  forward  biased,  a  net  imbalance 
in  the  currents  results  that  will  couple  to  the  output  guide  with  a 
reference  phase  of  zero  degrees.  If  now  the  diode  biases  are  reversed, 
an  imbalance  of  the  currents  in  the  opposite  direction  will  result  in 
the  same  magnitude  of  coupled  signal  but  with  the  opposite  polarity. 

The  plus  or  minus  quadrature  signal  can  be  excited  as  above  with  a 
second  slot  spaced  one-quarter  guide  wavelength  from  the  first  slot,  as 
illustrated. 

Figure  C-2  illustrates  the  eight  possible  quantized  states  for  the 
variphase  exciter.  Four  phase  states  are  realized  by  +1  and  +p.  An 
additional  four  phase  states  are  realized  with  the  combined  states  of 
the  device.  A  total  of  eight  phase  states  are  therefore  achievable  with 
only  four  diodes.  However,  the  coupling  in  the  diagonal  states  is  nomin¬ 
ally  3  dB  stronger  than  in  the  cardinal  state. 

The  absolute  value  of  coupling  to  the  output  element  is  controlled 
by  the  slot  depth.  Therefore,  only  one  design  of  switching  network 
comprising  the  printed-circuit  exciter  is  required  for  a  wide  range  of 
antenna  designs. 

Table  c-I  is  a  summary  of  the  projected  performance  of  a  linear 
array  utilizing  the  variphase  exciter.  The  array  has  110  elements  and 
a  nominal  beamwidth  of  one  degree.  The  phase-randomization  function  is 
quadratic  with  an  edge  phase  value  of  360  degrees,  or  one  wavelength  of 
curvature.  The  overall  beam  position  accuracy  and  sidelobe  performance 
are  similar  to  a  conventional  array  with  three-bit  phase  shifters. 


24 


Fig.  C-l  —  Vari phase  exciter  (narrow  wall  edge  slots). 


IN-PHASE  COMPO - 


Fig.  C-2  —  Quantized  states  of  variphase  exciter. 
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Table  C-I 


Projected  performance  of 

a  110-element  linear 

array 

Nominal 

Scan  Angle 
(Degrees) 

Beam  Position 
Error  (Degrees) 

Peak  Sidelobe 

Level  (dB) 

RMS  Sidelobe 
Level  (dB) 

0 

0.00 

-23 

-31 

1 

.01 

-22 

-31 

2 

0 

o 

o 

-23 

-31 

3 

-.01 

-22 

-31 

4 

O-l 

-24 

-31 

5 

-oOl 

-21 

-30 

6 

.00 

-22 

-31 

7 

.00 

-24 

-31 

8 

-.01 

-21 

-31 

9 

1 

• 

o 

H4 

-23 

-32 

10 

-.01 

-23 

-31 

15 

.01 

-24 

-31 

20 

.00 

-24 

-31 

25 

.00 

-22 

-31 

30 

-.02 

-21 

-30 

35 

.01 

-22 

-30 

40 

.03 

-22 

-30 

45 

-.01 

-21 

-30 
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2.  THE  COMPACT  ANTENNA 


Figure  C-3  shows  chat  Che  ideal  COMPACT  element  pattern  is  obtained 
by  an  element  that  produces  an  aperture  excitation  that  is  in  the  form 
of  sin  x/x.  The  width  of  the  "main-lobe"  of  the  distribution  in  wave¬ 
lengths  and  the  sharpness  of  the  pattern  cut-off  is  proportional  to  the 
antenna  aperture  size,  typically  on  the  order  of  a  few  degrees  wide. 

Figure  C-4  shows  a  block  diagram  of  a  typical  COMPACT  antenna.  It 
consists  of  a  power  divider  network,  one  phase  shifter  per  effective 
element,  a  coupling  network  and  individual  radiators.  The  power  divider 
network  is  used  to  excite  the  element  terminals  with  a  tapered  distribu¬ 
tion  that  yields  an  acceptable  sidelobe  level  for  the  scanning  beam 
pattern.  The  phase  shifters  shown  at  the  element  terminals  are  exercised 
to  scan  the  narrow  beam.  Inserted  between  the  element  terminals  and 
the  radiators  in  the  aperture  is  a  coupling  network.  This  network  is 
a  key  component  whose  function  is  to  provide  the  sin  x/x  element  termin¬ 
als,  The  network  must  provide  translational  symmetry;  that  is,  each 
element  terminal  excites  the  same  sin  x/x  distribution  in  the  aperture 
only  translated  in  the  aperture  by  the  separation  of  the  element  termin¬ 
als. 


Before  describing  the  actual  COMPACT  coupling  network,  the  basic 
method  of  generating  sin  x/x  with  translational  symmetry  is  shown  in 
Figure  C-5 .  The  network  synthesis  starts  with  a  horizontal  coupling 
line  and  vertical  element  lines  connected  to  it  by  identical  directional 
couplers,  A  portion  of  a  signal  input  at  one  of  the  element  lines  is 
directly  transmitted  to  the  output  of  that  element  line  and  the  remain¬ 
ing  portion  then  excites  all  element  lines  to  the  right  of  the  input  by 
means  of  the  directional  couplers.  This  repeated  coupling  process 
results  in  an  amplitude  distribution  which  monotonically  decreases  to 
the  right  of  the  element  terminal  input.  The  phasing  along  the  coupling 
line  between  directional  couplers  is  tt  radians  (180°  of  phase  shift). 

This  va'ue  of  phase  is  chosen  to  produce  the  correct  in-phase  polarity 
of  the  output  distribution  mainlobe  and  the  alternating-phase  polarity 
of  the  sidelobes,  as  in  the  sin  x/x  function.  As  shown  in  the  figure, 
the  inherent  property  of  a  directional  coupler  is  that  the  coupled  out¬ 
put  is  in  quadrature  ( j )  (90°  of  phase  shift)  with  the  directly  trans¬ 
mitted  output.  By  tracking  the  signal  input  through  this  basic  network, 
it  can  be  seen  that  the  proper  polarities  are  achieved.  This  network, 
however,  generates  only  one  half  of  the  desired  aperture  excitation,  and 
the  sampling  points  of  sin  x/x  are  rather  widely  spaced.  To  generate 
the  complete  sin  x/x  function  and  with  more  sample  points,  two  identical 
coupling  lines  are  interleaved  and  combined  at  the  bottom  of  the  network 
to  form  the  element  terminal.  The  directional  couplers  in  the  top 
coupling  line  excite  one  side  of  the  aperture.  The  directional  couplers 
in  the  bottom  coupling  line  excite  the  opposite  side  of  the  aperture. 

The  couplers  are  identical  in  both  coupling  lines  only  oriented  in 
opposite  directions.  By  tracking  signals  through  this  network,  it  can  be 
seen  that  an  element  terminal  excites  a  sin  x/x  type  distribution  in 
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(b)  ELEMENT  APERTURE  EXCITATION 

Fig.  C-3  —  Ideal  COMPACT  element  pattern  and  aperture  excitation. 
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Fig.  C-4  —  Block  diagram  of  a  COMPACT  antenna. 
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Fig.  C-o  —  Principle  of  COMPACT  coupling  network. 


the  aperture  with  four  samples  within  the  main  lobe,  one  sample  at  each 
of  the  zeros,  and  one  sample  at  the  peak  of  each  sidelobe  This  net¬ 
work  is  representative  of  the  coupling  network  required  in  COMPACT. 

The  network  of  Figure  C-6  has  the  required  translational  symmetry 
to  form  an  array  antenna.  Each  element  terminal  excites  the  same  aper¬ 
ture  distribution  only  displaced  by  the  element  spacing  because  the  net¬ 
work  is  identical  as  seen  from  each  element  terminal. 

The  implementation  for  coupling  networks  is  shown  in  Figure  C-7  . 

The  network  previously  described  is  used  with  a  slight  modification. 

This  modification  is  the  location  of  resistive  components  between  each 
directional  coupler  along  the  coupling  lines.  Each  attenuator  is  iden¬ 
tical,  as  are  the  couplers,  and  its  main  purpose  is  to  provide  one  addi¬ 
tional  degree  of  control  of  the  element  aperture  excitation.  This 
results  in  a  closer  approximation  to  the  desired  sin  x/x  and  closer 
realization  of  the  ideal  element  radiation  pattern  The  value  of  cou¬ 
pling  and  resistive  components  are  both  approximately  -3  dB  with  a  net¬ 
work  loss  of  about  -2  dB. 
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O  INTERLEAVE  TWO  OPPOSITE  COUPLING  LINES 
°  COMBINE  ELEMENT  LINES  IN  PAIRS  TO  FORM  TERMINAL 
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Fig.  C-6  -  Typical  COMPACT  coupling  network. 


APERTURE 


WAVE 

EXCITES 

RIGHT 

SIDE 


33 


MARGIN  FOR  FINITE  ANTENNA  BEAMWIDTH  AND  ELEMENT 
PATTERN  CUTOFF  WHEN  BEAM  SCANNED  TO  0max  AND 
®min 


Fig.  C-"  -  Implementation  of  the  COMPACT  coupling  network. 


APPENDIX  D 


Chronology  of  Events  and  Milestones  of  NRL 
Participation  in  MLS  under  Interagency  Agreement 

Date  Description 

1975 

May  Initiation  of  interagency  agreement* 

June  COMPACT  proprietary  proposal  by  Hazeltine. 

Meetings  with  Australians  and  French  in  Paris  and 
Sydney  * 

July  Proprietary  proposal  by  ITT  Gilfillan, 

August  DRB  meeting  at  Bendix,  Towson,  MD. 

September  DRB  meetings  at  TI,  Dallas  and  Bendix,  Ft.  Lauderdale. 
October  Preparation  of  contributions  to  U.S.  ICAO  submission. 
November  Visit  by  Russian  delegation  to  U.S. 

DRB  meeting  at  Bendix,  Towson,  MD. 

December  Meeting  with  Australians  in  Sydney. 

Contract  award  to  ITT  Gilfillan. 

1976 

February 

March 
April 
May 


June 
July 

Augus  t 
September 
November 

1977 

January  Review  of  NASA  Basic  Wide  proposals. 

March  AWOP  meeting  in  Montreal. 

April  Contract  award  to  Hazeltine  for  testing  of  COMPACT 

antenna  at  NAFEC. 

Performance  simulation  begun. 

June  Proposal  of  Hadamard  monitoring  algorithm. 


Contract  award  to  Hazeltine. 

(WGA  meeting  in  Germany  -  no  NRL  participation.) 

DRB  meeting  at  Bendix,  Towson,  MD. 

Completion  of  thinned  array  analysis. 

Phased  array  background  paper  completed. 

ITT  Gilfillan  contract  complete. 

WGA  meeting  in  Washington. 

Thinned-array  test  program  at  NAFEC. 

(WGA  meeting  at  The  Hague  -  no  NRL  participation.) 
Hazeltine  contract  complete. 

Phased  array  monitoring  simulation  begun. 

Dispute  over  delivery  of  Hazeltine  COMPACT  antenna. 
WGA  meeting  in  London. 
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1977 


June  Completion  of  phased  array  monitoring  report. 

October  Phase  shifter  measurements  begun. 

1978 


May 

July 

September 

November 


Proposal  of  3D  bootlace  lens  techniques. 
Proposal  of  optimum  beam  steering  technique. 
Visit  to  DOT  and  CSIRO  in  Australia. 
Presentation  of  '360-deg  antenna  design  study  to 
visiting  Russian  delegation. 


1979 


May  Visit  to  Russia. 

July  Contract  award  to  Myers  Associates. 
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APPENDIX  C. 

PERFORMANCE  SIMULATION  OF  MLS  PHASED  ARRAYS 


1.  INTRODUCTION 

The  Federal  Aviation  Administration's  (FAA)  Microwave  Landing 
System  (MLS)  uses  a  Time  Reference  Scanning  Beam  (TRSB)  system.  This 
TRSB  system  is  an  air-derived  data  system  which  operates  at  C-band 
(5030-5090  MHz)  to  provide  precision  landing  guidance  information  on  an 
aircraft's  azimuth  and  elevation  positions.  The  angular  position  of  the 
aircraft  is  determined  by  receiving  separate  azimuth  and  elevation  fan 
beams  generated  from  two  linear  phased  arrays  (in  the  cases  treated 
here)  which  scan  "TO"  and  "FRO"  across  the  coverage  sector  in  both  azi¬ 
muth  and  elevation  directions.  The  airborne  subsystem  measures  the 
elapsed  time  between  a  given  "TO"  and  "FRO"  scan.  This  time  difference 
is  directly  proportional  to  the  aircraft's  angular  position  (see  Fig. 
E-l).  In  order  to  determine  the  aircraft's  precise  angular  location, 
it  is  of  utmost  importance  to  measure  accurately  the  time  the  main 
beam  passes  the  aircraft.  The  stringent  MLS  operational  requirements 
set  this  position  measurement  at  an  accuracy  of  +.02  degree  with  a  one 
degree  beamwidth  fl].  The  accuracy  of  this  measurement  depends  on  both 
the  antenna  beam  shape  and  the  method  used  to  detect  the  peak  of  the 
beam  [2].  The  dynamic  antenna  pattern  is  affected  by  both  the  multipath 
reflections  and  the  array  radiation  pattern.  The  multipath  effect  in 
turn  depends  on  the  environment  of  the  antenna  site,  the  array  design 
and  the  beam  peak  detection  algorithm. 

A  computer  simulation  has  been  developed  which  represents  the  per¬ 
formance  of  the  various  phased  array  antenna  systems  which  are  being 
considered  for  use  in  the  Microwave  Landing  System.  The  purpose  of  this 
simulation  is  to  assess  the  angular  accuracy  of  MLS  in  terms  of  its 
design  parameters,  tolerance  effects,  component  and  subsystem  failures, 
multipath,  and  angle  detection  algorithm. 

In  this  way  it  is  hoped  that  a  basis  can  be  established  for  speci¬ 
fying  parameter  values  and  tolerances  for  the  ground  scanning  antennas 
and  for  determining  its  failure  modes  and  vulnerability.  Important  appli¬ 
cations  of  this  process  are  the  phase  and  amplitude  tolerances  on  the 
antenna  array,  which  are  in  turn  related  to  the  specifications  on  the 
phase  shifters  and  distribution  networks. 

Once  the  design  parameters  and  tolerances  and  the  failure  modes 
are  specified,  another  simulation  which  involves  the  fault  location 
monitoring  system  can  be  used  to  determine  whether  an  out-of-limits 
accuracy  condition  can  be  detected  in  a  timely  fashion. 
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In  the  following  sections  of  this  appendix,  the  simulation  of  the 
antenna  array  and  its  feed  system  is  discussed.  Then  the  computation 
of  the  radiation  patterns  is  described.  The  measurement  of  angle  oy 
the  airborne  receiver  is  analyzed.  The  effects  of  receiver  location 
and  multipath  are  treated.  Finally,  the  various  simulation  output  for¬ 
mats  are  presented,  and  a  flow  diagram  for  the  entire  simulation-is 
given  with  instructions  for  using  it. 

2.  SIMULATIONS  OF  ARRAYS 

This  program  simulates  the  performance  of  linear  phased  arrays 
used  in  the  TRSB  system.  Two  types  of  arrays  are  simulated.  The  first 
type  is  a  conventional  phased  array  in  which  each  radiating  element  is 
equipped  with  a  beam  steering  phase  shifter.  The  second  type  is  a 
COMPACT  array.  Several  antenna  radiating  elements  in  this  array  may 
share  a  common  phase  shifter.  The  operational  principles  of  this  array 
are  included  in  Appendix  C.  Array  data  such  as  element  spacing,,  number 
of  phase  shifters,  number  of  radiating  elements  per  phase  shifter  and 
the  phase  shifter's  number  of  bits  are  all  variable  inputs  to  the  sim¬ 
ulation  program.  The  maximum  number  of  radiating  elements  in  this 
simulation  program  is  limited  to  120  while  the  number  of  phase  shifter 
bits  must  be  no  more  than  6. 

It  is  assumed  that  both  errors  and  failures  exist  in  the  array. 

In  the  following,  simulations  of  errors  and  failures  are  described. 

(a)  Errors: 

All  components  in  the  array,  which  includes  the  feed  network,  rf 
power  source  and  phase  shifters,  contribute  both  amplitude  and  phase 
errors.  In  the  simulation  program,  errors  from  all  sources  are  lumped 
together  and  assigned  to  each  phase  shifter.  The  distribution  of  these 
errors  is  assumed  to  be  Gaussian  in  view  of  the  fact  that  they  are  con¬ 
tributed  from  many  independent  sources. 

Amplitude  and  phase  errors  are  assumed.  Amplitude  of  radiating 
power  from  each  phase  shifter  is  assumed  normally  to  be  unity.  Due  to 
feed  network  variation  and  phase  shifter  attenuation,  amplitude  of  rf 
power  varies  from  element  to  element.  This  variation  is  assumed  to  be 
Gaussian  centered  at  unity.  The  standard  deviation,  in  percent,  of 
this  variation  must  be  specified.  As  for  phase  error,  two  types  are 
assumed.  The  first  one  is  an  insertion  phase  error  which  is  assumed  to 
be  associated  with  each  phase  shifter,  and  the  second  one  is  a  phase 
shifter  error  which  is  associated  with  each  phase  state  of  each  phase 
shifter.  The  former  varies  from  element  to  element  while  the  latter 
varies  from  phase  state  to  phase  state  of  each  phase  shifter.  Both  of 
these  errors  are  assumed  to  be  Gaussian  centered  at  a  normal  value. 
Standard  deviations  of  these  errors  are  inputs  to  the  program.  Besides 
these,  a  constant  insertion  phase  can  be  added  to  each  phase  shifter 
across  the  array. 


For  a  COMPACT  array,  additional  phase  and  amplitude  errors  are 
assigned  to  the  coupling  network.  Again  these  errors  are  assumed  to  be 
Gaussian  centered  at  the  normal  value  with  a  prescribed  standard  devi¬ 
ation.  For  amplitude  error,  its  standard  deviation  is  in  percent  while 
phase  errors  are  stated  in  degrees.  These  standard  deviations  are  input 
to  the  simulation  program.  A  subprogram  which  generates  Gaussian  ran¬ 
dom  numbers  with  zero  mean  and  unit  standard  deviation  is  included  in 
the  simulation  program.  Outputs  from  this  subprogram,  multiplied  by 
the  appropriate  standard  deviation  value  are  used  as  the  simulated  errors . 
At  the  beginning  of  each  sample  of  simulation,  the  amplitude  error  and 
the  insertion  phase  error  associated  with  each  phase  shifter  and  the 
phase  error  of  each  phase  state  are  generated  and  are  stored  in  accor¬ 
dance  with  the  location  and  phase  state  of  each  phase  shifter.  For  a 
COMPACT  array,  errors  of  amplitude  and  phase  in  the  coupling  network  are 
stored  in  accordance  with  the  output  port  of  the  network.  These  recorded 
values  are  used  for  computation  of  the  radiating  energy  level  received 
by  the  MLS  receiver. 

(b)  Failures: 

The  number  of  failures  is  an  input  parameter  to  the  simulation 
program.  These  failures  can  occur  on  any  phase  shifter  element  and  at 
any  bit  of  that  phase  shifter  element.  Three  modes  of  failures  are 
assumed.  The  first  one  is  an  open  circuit  failure.  This  may  be  taken 
into  account  for  any  point  along  the  feed  network  from  rf  source  to  the 
radiating  element.  A  prescribed  probability,  an  input  to  the  simula¬ 
tion  program,  is  assigned  to  this  open-circuit  failure.  The  second  and 
third  types  of  failures  are  failures  of  the  phase  shifter  element  itself . 

The  second  type  is  a  failure  such  that  a  phase  shifter  may  stick  to 
a  certain  phase  state  and  fail  to  respond  to  the  beam  steering  command. 
Probability  of  occurrence  of  this  type  of  failure  is  assumed  to  be  known. 
The  third  type  is  a  failure  of  one  of  the  phase  shifter  bits.  This  may 
be  a  stuck-to-zero  or  a  stuck-to-one  failure.  The  probability  of  any 
phase  state  in  each  phase  shifter  to  develop  an  error  of  the  second 
type  is  assumed  to  be  equal.  Similarly,  all  bits  in  a  phase  shifter 
have  equal  probability  to  develop  either  a  stuck-to-zero  or  a  stuck-to- 
one  failure. 

A  subprogram  which  generates  uniformly  distributed  random  numbers 
is  included  in  the  simulation  program.  By  use  of  the  random  numbers 
generated  from  this  subprogram,  locations  and  types  of  failures  are 
determined.  The  procedure  for  performing  such  a  task  in  the  simulation 
program  is  as  follows. 

For  each  failure  specified,  a  random  number  is  drawn  first.  Since 
this  random  number  has  a  uniform  distribution  from  0  to  1,  multiplying 
this  random  number  by  the  number  of  array  elements  and  truncating  the 
fraction  portion  can  determine  the  location  of  a  randomly  failed  ele¬ 
ment.  By  the  same  token,  a  second  random  number  can  be  used  to  determine 
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Che  type  of  failure.  The  random  number  which  has  a  range  from  0  to  1 
is  divided  into  three  regions.  Each  region  represents  a  type  of  failure. 
For  example,  if  we  specify  .1  for  open  circuit  failure  and  .2  for  stuck- 
phase-state  failure,  then  random  numbers  lying  between  zero  and  .1  repre¬ 
sents  open  circuit  failure,  .1  to  .2  represents  stuck-phase-state  fail¬ 
ure,  and  greater  chan  .2  represents  the  case  of  bit  failure.  A  third 
random  number  is  then  used  to  determine  which  bit  or  which  phase  state 
is  failing. 

Additional  failures  may  occur  in  the  coupling  network  of  a  COMPACT 
array.  These  failures  are  assumed  to  be  open  circuit.  They  may  occur 
at  any  input  port,  output  port  or  coupling  line  with  equal  probability. 
Locations  of  these  failures  can  be  determined  by  the  uniform  random 
number  generator.  The  number  of  failures  in  the  COMPACT  array  must  also 
be  specified. 

Distributions  of  these  errors  and  failures  are  generated  at  the 
beginning  of  each  simulation  sample.  In  order  to  achieve  some  statis¬ 
tical  distribution  of  the  effects  of  these  errors  and  failures  on  the 
performances  of  the  array,  a  number  of  simulation  samples  usually  are 
performed.  However,  among  all  these  samples,  the  number  of  failures, 
the  probability  of  certain  types  of  failures  and  the  standard  devia¬ 
tions  of  errors  are  maintained  at  a  constant  value.  The  number  of  simu¬ 
lation  samples  is  an  input  variable,  and  is  specified  at  the  beginning 
of  the  program. 

3.  ARRAY  PATTERN  COMPUTATION 

For  each  simulation  sample,  a  number  of  receiver  locations  are 
allowed.  At  each  receiver  location,  the  direct  phase  path  and  the 
reflection  phase  path  (if  multipath  is  assumed)  from  each  radiating  ele¬ 
ment  to  the  receiver  are  computed  and  stored.  For  most  cases,  except 
for  the  monitor  output,  the  receiver  is  assumed  to  be  located  in  the  far 
zone  region.  Hence,  the  phase  path  is  a  function  only  of  the  receiver ' s 
location  angle.  This  phase  can  be  represented  as 

P(9,n)  =  exp  [j2TT'nd/X  sin  0], 

where  9  is  the  receiver  position  angle  with  respect  to  the  normal  to  the 
plane  of  the  array,  n  is  the  index  of  the  array  element  and  d  is  the 
element  spacing.  The  reflection  phase  path  depends  on  the  nature  of  the 
reflection  surface.  This  will  be  discussed  in  a  later  section.  For  the 
case  of  the  monitor  output,  the  monitor  range  is  small,  and  its  phase 
path  is  calculated  based  on  the  actual  distance  from  monitor  receiver  to 
the  array,  which  can  be  represented  as 

D2  =  (H  -H  +nd) 2  +  R2 
n  2  l 

P(9,n)  *  exp  [j  2ttD  /X], 
n 
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where  and  R  are  the  receiver  height  and  range,  respectively,  and  H1 
is  the  antenna  height  at  the  center  of  the  array. 

The  array  beam  scans  from  one  side  of  the  receiver  to  the  other. 
This  constitutes  a  "TO"  scan,  and  reversed  scan  constitutes  a  "FRO" 
scan.  The  scan  algorithm  or  beam  steering  algorithm  is  predetermined 
and  is  generated  in  a  subprogram  that  will  be  described  in  a  later 
section.  In  order  to  reduce  simulation  computation,  received  signals 
are  computed  only  in  a  small  sector  which  is  approximately  equal  to  one 
beamwidth  centered  at  the  receiver's  angle  position.  At  each  beam 
steering  step,  the  required  phase  settings  are  first  computed  in  accor¬ 
dance  with  the  beam  steering  algorithm.  The  previously  stored  phase 
state  settings,  which  include  both  amplitude  and  phase  errors  and  fail¬ 
ures,  are  then  called  and  used  to  compute  the  field  radiating  from  each 
element  in  the  direction  of  the  receiver.  This  radiating  field  is  then 
translated  through  both  the  direct  path  and  the  reflection  phase  path 
from  the  element  to  the  receiver.  The  total  field  reaching  the  receiver 
at  this  beam  position  is  the  sum  of  the  field  components  from  the  radi¬ 
ating  elements.  This  total  radiated  field  at  each  beam  steering  step 
in  both  "TO"  and  "FRO"  directions  is  stored  for  angle  processing  by  the 
receiver. 

A  Taylor  weighted  illumination  subprogram  is  included  in  the  simu¬ 
lation  program.  One  may  specify  the  desired  sidelobe  level  and  the 
design  parameter  n.  The  subprogram  generates  the  required  array  illum¬ 
ination  function. 

For  a  COMPACT  array,  besides  tapered  illumination,  the  program  also 
provides  variations  of  the  coupling  coefficient  and  attenuation  in  the 
"COMPACT"  coupling  network,  tilting  of  the  element  pattern  to  a  desired 
angle  and  the  effect  of  changing  the  operating  frequency.  Dummy  ele¬ 
ments  can  be  added  at  each  end  of  the  array. 

4.  ANGLE  RESOLUTION 

The  stored  dynamic  pattern  data  is  first  smoothed  by  a  two-pole 
digital  Butterworth  filter.  Details  of  the  filter  design  are  given  in 
Appendix  F.  These  smoothed  data  are  then  used  for  angle  resolution. 

There  are  several  methods  to  determine  the  receiver  location  angle. 
The  one  used  in  this  simulation  program  is  that  proposed  by  Bendix  in 
which  Che  two  half  power  points  of  the  scanning  dynamic  pattern  are 
first  located  [3].  The  center  of  the  beam  is  taken  as  the  average  of 
these  two  points.  In  the  actual  system,  the  array  beam  scans  continu¬ 
ously  in  the  coverage  sector  from  -0j  to  +0i  (see  Fig.  E-l).  The  air¬ 
borne  receiver  senses  the  radiated  energy,  smooths  it,  determines  the 
3-dB  points  and  records  the  time  of  the  middle  point,  as  it  passes  the 
aircraft.  The  time  difference  between  a  "TO"  and  "FRO"  scan  is  used 
to  determine  the  aircraft's  angle  location.  However,  in  a  simulation, 
the  amount  of  computation  required  to  scan  over  the  entire  coverage 
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sector  is  excessive.  Therefore,  in  the  simulation  program  the  dynamic 
pattern  is  simulated  only  in  the  main  beam  region.  Since  the  beam  scans 
in  discrete  steps,  each  beam  position  in  the  coverage  section  is  indexed 
in  the  simulation  program  in  accordance  with  the  beam  steering  algorithm. 
It  is  further  assumed  that  for  both  the  "TO"  and  "FRO"  scans,  the  beam 
starts  precisely  at  the  same  beam  positions.  In  the  simulation  program, 
all  angle  positions  are  referenced  to  the  beam  indices  in  order  to 
achieve  better  accuracy  and  avoid  computer  truncation  error.  Two  prob¬ 
lems  have  to  be  addressed.  The  first  is  that  due  to  multipath  effects, 
the  beam  shape  may  be  deformed  and  the  main  beam  may  be  widened.  Start¬ 
ing  the  beam  at  the  half-beamwidth  point  of  an  ideal  pattern  may  not 
include  the  actual  3-dB  points.  To  overcome  this  difficulty,  an  adaptive 
algorithm  is  built  into  the  simulation  program.  The  dynamic  pattern  is 
first  computed  taking  into  account  the  effect  of  multipath,  errors  and 
failures.  The  3-dB  points  are  then  searched.  If  no  3-dB  point  can  be 
found  at  one  side  of  the  dynamic  pattern,  that  side  of  the  pattern  is 
extended  and  recomputed.  This  extension  will  continue  until  either  a 
3-dB  point  is  reached  or  the  pattern  level  is  not  reducing  any  more. 
Another  problem  is  the  transient  response  of  the  digital  filter.  To 
avoid  this  problem  the  computed  dynamic  pattern  must  have  enough  lead 
time  to  avoid  the  transient  region.  This  difficulty  is  avoided  by  first 
examining  the  delay  of  the  impulse  response  of  the  digital  filter.  In 
the  subsequent  pattern  computation,  the  range  of  the  pattern  is  extended 
to  cover  this  transient  delay  region.  By  use  of  both  the  "TO"  and  "FRO" 
patterns  the  effect  of  delay  due  to  filtering  is  removed. 

5.  RECEIVER  LOCATION 

For  each  simulation  sample,  a  number  of  different  receiver  loca¬ 
tions  may  be  simulated.  However,  in  all  these  simulations,  the  phase 
shifter  and  feed  network  errors  and  the  failures  remain  the  same.  The 
number  of  different  receiver  locations  to  be  simulated  is  an  input 
parameter.  There  are  five  ways  to  specify  these  receiver  locations. 

a.  The  receiver's  locations  are  specified  in  height  and  range. 

These  heights  and  ranges  are  read-in  to  the  program.  The  antenna  height 
at  array  center  must  also  be  specified.  The  location  of  the  monitor 
receiver  can  be  entered  in  this  way. 

b.  Both  receiver  height  and  range  are  generated  randomly  with  a 
uniform  probability  within  a  given  range.  Antenna  height  must  be 
specified. 

c.  In  most  cases,  the  receiver  is  in  the  far  field  of  the  antenna, 
and  only  the  angle  position  of  the  receiver  is  of  interest.  In  these 
cases  the  receiver  angle  locations  may  be  read  in. 

d.  Given  the  angle  ranges,  the  receiver  angle  locations  may  be  gen¬ 
erated  in  random  with  uniform  probability  within  a  specified  angle  range. 
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e.  For  pole  tests,  the  receiver  is  located  at  a  fixed  range  and  the 
height  of  the  receiver  is  increased  by  a  fixed  amount  until  a  certain 
height  is  reached  or  a  given  number  of  pole  test  positions  are  evaluated. 
For  this  case,  the  antenna  height  must  also  be  specified. 

6.  MULTIPATH  EFFECT 

The  multipath  effect  is  calculated  for  the  elevation  antenna  on  the 
assumption  that  the  ground  is  flat.  The  ground  reflection  angle  (see 
Fig.  E-2a)  is  then 


where  Hi  and  H2  are  the  antenna  height  and  receiver  height,  respectively, 
and  R  is  the  range.  The  phase  of  the  reflected  energy  is 

9p  =  4ttHiH2/RA+<|> 

Where  A  is  the  wavelength  and  <p  is  the  ground  reflection  phase  which 
along  with  the  ground  reflection  coefficient  must  be  specified.  For 
the  case  when  receiver  location  is  specified  only  by  elevation  angle, 
the  ground  reflection  angle  ip  -  -6,  where  9  is  elevation  angle.  The 
ground  reflection  phase  is  then  chosen  randomly  in.  the  range  from  0  to 
2tt. 


7.  BEAM  STEERING  ALGORITHM 

Four  beam  steering  algorithms  are  incorporated  into  the  simulation 
program.  However,  any  other  steering  algorithm  can  be  adopted  easily 
by  changing  a  subprogram  which  is  used  to  generate  these  steering 
algorithms.  The  four  beam  steering  algorithms  are  described  in  the 
following: 

a.  Bendix  algorithm  -  The  algorithm  was  developed  by  Bendix  for  a 
conventional  array.  In  this  algorithm,  the  array  phase  settings  are 
computed  at  larger  coarse  steps.  However,  from  one  coarse  step  to  the 
next  coarse  step,  the  phase  shifters  are  switched  one  pair  at  a  time. 
Thus  there  are  N/2  fine  steps  between  each  coarse  step,  where  N  is  the 
total  number  of  phase  shifters  (even)  in  the  array.  To  assure  uniform 
scan  from  one  coarse  step  to  the  next,  Bendix  has  developed  a  special 
sequence  to  achieve  this  goal.  That  sequence  is  used  in  this  simulation 
program. 


b .  Random  sequence  -  This  is  similar  to  that  of  the  Bendix  algo¬ 
rithm.  However,  the  switching  sequence  between  coarse  steps  is  random 
and  is  generated  in  accordance  with  a  random  number  sequence. 

c.  Hazeltine  sequence  -  This  sequence  wars' developed  by  Hazeltine  for 
use  in  their  COMPACT  array.  In  this  sequence  a  single  phase  shifter  at 
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a  time  is  switched.  The  interval  of  switching  time  is  2/7  usee.  For 
an  array  scan  rate  of  20,000  degrees  per  second,  this  gives  a  scan  step 
of  approximately  .0057  degree.  The  Hazeltine  sequence  steers  the  beam 
from  -1°  to  10  degrees. 

d.  Conventional  step  -  In  this  straightforward  procedure,  for  the 
required  steering  angle  ac  any  step,  the  necessary  phase  settings  at 
that  step  are  computed  and  switched  accordingly. 

8.  SIMULATION  OUTPUT 

Several  different  kinds  of  simulations  can  be  implemented  with  this 
program.  These  simulations  and  their  outputs  are  described  in  the  fol¬ 
lowing  paragraphs  and  an  example  of  each  is  presented. 

a.  Angle  error  statistic  -  The  main  purpose  of  this  implementation 
is  to  establish  a  statistical  correlation  between  angle  error  and  com¬ 
ponent  tolerance  (array  errors)  and  the  number  of  failures  in  the  array. 
The  output  of  this  simulation  is  a  plot  of  the  cumulative  probability  of 
angle  error,  given  the  number  of  failures  and  the  standard  deviations  of 
amplitude  and  phase  errors.  In  the  example  of  Fig.  E-3,  the  probability 
distributions  of  angle  errors  for  a  conventional  96-element  array  are 
presented.  The  array  elements  are  subject  to  a  Gaussian  amplitude  and 
phase  error  with  standard  deviations  of  .05  for  amplitude  and  5  degrees 
for  phase.  Three  curves  are  presented,  one  for  zero  failure,  one  for  6 
failures  and  the  third  for  a  12  failure  case.  For  example,  when  the 
array  has  no  failure,  the  probability  that  the  angle  error  is  less  than 
.02  degrees  is  95%,  while  for  a  6  failure  case,  this  becomes  91%  and  for 
12  failures  it  drops  to  about  80%. 

This  statistic  is  generated  by  performing  a  number  of  simulation 
samples.  For  each  sample,  a  number  of  receiver  locations  are  chosen. 

In  this  case,  these  receiver  locations  are  chosen  at  random  with  a 
coverage  section  of  -60°  to  +60°.  In  this  particular  type  of  simula¬ 
tion  the  number  of  error/ failure  distributions  is  30  while  the  number 
of  receiver  locations  is  40.  This  gives  a  statistical  sample  of  1200. 

The  curves  of  Fig.  E-3  are  generated  by  these  samples. 

b.  Pole  test  -  the  pole  test  is  used  in  the  field  to  measure  the 

system  performance.  In  this  test,  a  receiver  is  located  at  a  given 

distance  from  the  array  antenna,  and  a  series  of  measurements  is 

made  at  given  increments  of  height.  The  output  of  these  simulations  is 
shown  in  Fig.  E-4a.  This  plot  is  Tor  the  case  of  a  COMPACT  array  with 
24  phase  shifters.  No  failure  is  assigned  to  the  array.  Five  percent 
and  five  degree  standard  deviations  are  assigned  respectively  to  the 
amplitude  and  phase  distributions  of  the  array  errors.  The  receiver  is 

at  a  range  of  280.4  meters  (920  ft),  starting  at  an  initial  height  of 

.61  meter  (2  ft)  and  increasing  in  .61  meter  (2  ft)  steps  until  it 
reaches  20.7  meters  (68  ft).  At  each  step  16  simulations  are  performed. 
For  each  simulation,  the  array  failures  and  errors  remain  the  same;  how¬ 
ever,  each  phase  shifter  is  incremented  by  one  phase  state.  This  process 
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is  known  as  phase  cycling.  Angle  errors  for  each  simulation  at  each 
receiver  height  are  then  computed.  The  mean,  standard  deviation,  mini¬ 
mum  and  maximum  of  this  angle  error  are  then  plotted  for  each  receiver 
height.  It  can  be  seen  that  due  to  the  multipath  effect,  angle  errors 
are  very  large  at  lower  receiver  heights.  Furthermore,  because  the 
element  pattern  of  this  array  is  tilted  12  degrees  upwards,  a  bias  angle 
error  exists.  In  this  example,  the  increment  of  the  phase  shifter  is 
set  at  one  bit  (360/16  degrees  for  a  4  bit  phase  shifter).  However,  the 
program  has  provision  to  accommodate  any  odd  increment  of  insertion 
phase. 

Hazeltine  developed  a  beam  steering  algorithm  which  can  effectively 
compensate  this  bias  error.  •  One  of  the  examples  is  shown  in  Fig.  E-4b. 

In  this  figure,  the  pole  test  was  performed  in  the  elevation  range  from 
2.25  to  3.75  degrees  with  11  angle  positions  evenly  distributed  within 
this  range.  At  each  angle  position,  10  samples  are  computed.  Each 
sample  has  phase  cycling  similar  to  the  case  discussed  above.  Two  sets 
of  data  are  presented  side  by  side  at  each  angle  position.  The  circles 
are  the  mean  value  of  the  16  phase  cycling  results.  At  each  angle 
position,  there  are  10  such  samples.  To  the  right  of  these  samples  are 
the  mean,  standard  deviation,  minimum,  and  maximum  of  these  samples. 

The  global  mean  is  the  average  value  of  all  the  mean  values  at  each 
angle  position,  while  standard  deviation  is  the  mean  value  at  each  angle 
position  deviated  from  this  global  mean.  The  bias  angle  errors  are 
seen  to  be  removed  for  this  example. 

c.  Dynamic  pattern  -  In  order  to  examine  closely  the  effects  of 
multipath,  beam  steering  algorithms  and  the  digital  filter,  the  scanning 
dynamic  pattern  is  plotted.  This  pattern  is  the  radiated  energy  received 
by  the  receiver  at  a  fixed  position.  One  of  the  examples  is  shown  in 
Fig.  E-5.  This  is  a  dynamic  pattern  of  a  receiver  at  the  horizon  for  a 
conventional  array.  There  are  four  curves  on  this  plot.  Two  of  the 
curves  which  are  almost  coincident  are  the  dynamic  pattern  for  both  "TO" 
and  "FRO"  cases  before  filtering.  The  two  curves  with  cross  marks  are 
the  patterns  after  filtering.  Delays  due  to  filtering  for  "TO"  and  "FRO" 
cases  are  in  opposite  directions.  One  may  see  that  at  angles  which  are 
below  the  horizon  (below  receiver's  height),  the  pattern  shapes  are 
strongly  modified  by  the  ground  reflections.  Furthermore,  no  3-dB  point 
exists  at  the  lower  angle  side,  even  with  the  pattern  extended  further 

as  shown  in  this  figure. 

d.  Monitor  output  -  The  monitor  output  is  a  plot  of  the  dynamic 
pattern  for  a  monitor  receiver  situated  at  close  range.  Because  the 
receiver  is  located  at  a  closer  range,  the  pattern  is  not  computed  on 
the  assumption  that  the  receiver  is  at  infinity.  Instead,  the  actual 
distance  from  each  element  to  the  receiver  is  computed  and  the  actual 
phases  of  both  the  direct  radiation  and  the  reflection  from  the  ground 
surface  are  used  to  compute  the  dynamic  pattern.  An  example  for  a 
COMPACT  array  is  shown  in  Fig.  E-6.  In  this  plot  the  mulipath  effect 
is  ignored. 
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e.  Array  or  element  pattern  -  Sometimes  it  is  desirable  to  be 
able  to  examine  the  array  pattern  or  the  element  pattern  (for  COMPACT 
array).  This  pattern  may  be  used  to  examine  the  effects  of  failures 
and  errors  in  the  array.  It  also  can  be  used  to  examine  the  effect  of 
element  pattern  tilt,  frequency  shifting  and  other  effects  on  the  ele¬ 
ment  pattern  of  a  COMPACT  array.  For  a  COMPACT  array  one  must  specify 
the  particular  element  in  the  array  to  be  plotted  and  also  the  pattern 
tilt  angle,  frequency  shifts,  etc.  One  example  is  shown  in  Fig.  E-7 . 

This  is  a  plot  of  the  element  pattern  of  a  COMPACT  array.  The  element 
is  at  the  center  of  the  array  of  116  radiating  elements  and  24  phase 
shifters  with  8  dummy  radiating  elements  located  at  each  end  of  the 
array.  The  element  is  tilted  at  a  12°  angle.  The  operating  frequency 
is  at  the  center  band. 

f.  Phase  shifter  switching  histogram  -  For  a  given  steering  algo¬ 
rithm  it  is  of  interest  to  examine  the  number  of  phase  shifters 
switched  at  each  coarse  step.  This  statistic  may  be  useful  in  design¬ 
ing  steering  algorithms.  The  beam  steering  algorithm  in  this  example 
is  that  proposed  by  Bendix.  The  coarse  step  is  .2  degree.  Fig.  E-8 
shows  the  number  of  phase  shifters  switched  at  each  coarse  step.  It  is 
interesting  to  note  that  about  one  third  of  the  phase  shifters  are  not 
switched  during  any  given  coarse  step. 

g.  Digital  filter  response  -  A  digital  filter  is  used  to  smooth  the 
dynamic  pattern  data  before  they  are  used  for  angle  resolution.  The 
response  of  this  filter  can  be  plotted  as  shown  in  Fig.  E-9.  Both  fre¬ 
quency  and  impulse  responses  are  shown  in  this  figure. 

9 .  FLOW  DIAGRAM 

A  flow  diagram  of  this  simulation  program  is  shown  in  Fig.  E-10. 

The  simulation  program  starts  with  a  read-in  of  all  pertinent  simulation 
data.  A  list  of  such  data  is  shown  in  Table  E-l.  If  the  simulation  is 
to  plot  the  smoothing  filter  response,  it  will  jump  to  the  subroutine 
for  generating  such  filter,  plot  its  responses  and  terminate  the  pro¬ 
gram.  Otherwise,  the  program  will  first  initiate  the  random  number 
generator.  Before  computing  the  radiating  scan  pattern,  the  program 
generates  the  beam  steering  algorithm,  sets  up  the  receiver's  angle 
positions,  computes  the  phase  and  amplitude  of  the  multiple  path  reflec¬ 
tion  wave,  generates  smoothing  digital  filter  and  computes  the  phase 
paths  of  array  elements  to  each  receiver  position.  The  program  has  three 
layers  of  do-loops.  The  outer  loop  is  the  sample  loop  in  which  the 
failure  modes  and  locations  and  the  array  errors  are  determined  for  each 
simulation  sample.  The  second  loop  is  the  phase  cycling  loop  in  which 
a  constant  phase  is  added  to  each  array  element  for  phase  cycling  pur¬ 
poses.  The  third  loop  performs  the  computation  of  the  dynamic  array  pat¬ 
tern  at  each  receiver  position.  After  that,  the  angle  position  will  be 
determined  and  compared  with  the  actual  angle  position.  Results  are 
then  plotted.  A  pre-scan  at  each  receiver  angle  position  is  performed  to 
make  sure  that  the  computed  dynamic  pattern  includes  the  3-dB  points  and 
avoids  the  transient  effect  of  the  smoothing  filter. 


oooooooooooooooooooooooooooooo ooooooooooooooooooooooooooooooo 


Table  E-l  —  Data  Read-in  Procedure  for  Simulation 


♦♦♦♦••♦FIRST  DATA  CARO,  USE  15  F0RRAT********«*«****************** 

NL,  NUMBER  OF  PHASE  SHIFTERS 

NM,  NUMBER  OF  RADIATING  ELEMENTS  PER  PHASE  SHIFTER 
NM=1 ,  CONVENTIONAL  ARRAY 

NBIT.  NUMBER  OF  PHASE  SHIFTER  BITS 
NSAMP,  NO.  OF  SIMULATION  SAMPLES 

NSAG,  NUMBER  OF  DIFFERENT  RECEIVER  ANGLE  POSITIONS 
NF,  NUMBER  OF  PLOTS  OF  DIFFERENT  NUMBER  OF  FAILURES, OR  DIFFERENT 
ELEMENT  ERRORS 

NF  POSITIVE  FOR  DIFFERENT  ELEMENT  ERRORS 

NF  NEGATIVE  FOR  DIFFERENT  NUMBER  OF  FAILURES 

LTALR  EQUAL  1  USE  TAYLOR  WEIGHT 

LPKI ,  INITIAL  PHASE  STATE  FOR  POLE  TEST 

IF  LPKI  .GT.  2**NBIT ,  PLOT  OUT  EACH  INSERTION  PHASE  STATE 
JPRT  PRINT  ALL  DATA  AT  JTH  ANGLE  POSTION  AT  THE  FIRST  SAMPLE  FOR 
DIAGNOSIS  PURPOSE 
LLCT=0 ,  NO  LABEL 

LREF=0  NO  REFLECTION  IS  CONSIDERED 

LPOLE=0  ,  NO  POLE  TEST 

LPOLE  GT  0,  PERFORM  POLE  TEST, NUMBER  OF  RUN  IN  EACH  TEST  EQUALS 
TO  LPOLE,  SAME  PHASE  AODEO  TO  AND  FRO  SCAN 
LPOLE  LT  0,  DIFFERENT  PHASE  IS  ADOED  AT  TO  AND  FRO  SCAN 

LPLT=0  SUPRRESS  ALL  PLOTS 

LPLT=1 ,  PLOT  ANGLE  ERROR  STATISTIC  OR  POLE  TEST 
LPLT=1 1 ,PLOT  ANGLE  ERROR  STATISTIC  WITH  PHASE  CYCLING, 

NO.  OF  OF  CYCLINGS  EQUAL  TO  LEAP 

LPLT=21 ,  PLOT  ONLY  MEAN  VALUE  OF  THE  PHASE  CYCLING 
LPLT=31 «  DO  POLE  TEST  FOR  MANY  SAMPLES  AND  PLT  ONLY  THE  MEAN  OF 
PHASE  CYCLING 

LPLT=41,  00  POLE  TEST  AT  ONE  SAMPLE  AND  PLOT  PHASE  CYCLING 
RESULTS 

LPLT=5 1 ,  PLOT  INDIVIDUAL  SAMPLE  FOR  SPACIAL  VARIATIONS 
LPLT=61 ,  00  3 OTH  LPLT=51,  AND  41 
LPLT=2,  PLOT  OYNAMIC  PATTERN 

LPLT=3,  PLOT  ARRAY  PATTERN  OR  ELEMENT  PATTERN 

LPL  T=4 ,  PLOT  PHASE  SHIFTER  SWITCHING  HISTOGRAM 
LPLT=  5 ,  PLOT  OUTPUT  OF  MONITOR 

LPLT=6,  PLOT  FILTER  FREQUENCY  AND  IMPULSE  RESPONSES 
INIT,  RANDOM  NUMBER  INITIALIZATION  SEED 

FIRST  DIGIT  OF  INIT  IS  USED  TO  SET  INITIAL  RANDOM  SEQUENCE 

2nd  Digit  multiplied  by  sample  index  sets  the  sample  sequence 

LEAP=0,NO  PATTERN  PLCT,LEAP=1,  PLOT  ELEMENT  PATTERN,  LEAP=2,  PLOT 
THE  FULL  ARRAY  PATTERN 

WHEN  L PLT =2 1,31, 41, 51, 61,  LEAP  IS  USED  AS  TH  INDEX  OF  THE 
REFERENCE  TO  COMPUTE  RMS  DEVIATION 
LSEQ,  BEAM  STEERING  SEQUENCE  CONTROL 
LSEQ=0 ,  NO  SEQUENCE 
LSEQ=1 »  USE  NAFEC  SEQUENCE 
LSEQ=2,  USE  RANDOM  SEQUENCE 
LSEQ=3 ,  USE  HAZETINE'S  SEQUENCE 

♦♦♦♦♦♦♦SECOND  DATA  CARD  USE  F10.6  FORNAT*******«***»*******»*«***» 
FAC,  COURSE  INCREMENT  CIN  DEGREES) 

CPA  ARRAY  TILTING  ANGLE 

REF  REFLECTION  COFFICIENT 

DX,  RADIATING  ELEMENT  SPACING 

RATIO  SCAN  PATTERN  SEARCH  CUT-OFF  POINT, PER  CENT  OF  MAX.  POWER 
OMEGA,  FILTER  CUT-OFF  FREQUENCY 
PTSR,  PLOT  SCALE  RANGE 

TAOLLT  SCAN  ANGLE  LIMITATION  BELOW  HORIZON 
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Table  E-l  —  (continued) 

......THIRD  OATA  CARO  USc  F10.6  FCRM4T,  SKIP  IT  IF  NOT  COMPACT  ARRAY 

TL  AG  ,  COMPACT  ARRAY  ELEMENT  PATTER#  TILTING  ANGLE 
-RAT,  PER  OF  NT  OF  FREQUENCY  DEVIATION  FROM  CENTER  FREQUENCY 
-A,  THE  FACTOR  THAT  THREE  D6  COUPLER  DEVIATES  FROM  UNIT  WHEN 
FREQUENCY  IS  CHANGED 

ENL*  THE  ELEMENT  THAT  THE  ELEMENT  PATTERN  TO  BE  PLOTTED 
XR,  YR  COMPENSATION  COE FF.  FOR  COMPACT  ARRAY 
-LV3C.  NUMBER  OF  DUMMY  ELEMENTS  IN  A  COMPACT  ARRAY 

............FOURTH  DATA  CARD  USE  F1C.6  FORMAT************************** 

Hit  ANTENNA  HEIGHT 

H?  TARGET  HEIGHT  OR  RANGE  OF  TARGET  HEIGHTS 
HE,  FOR  POLE  TEST  INCREMENT  OF  RECEIVER  HEIGHT 
cOR  MONITOR  ,  LIMIT  OF  LOWER  ANGLE 
DR  RANGE  OF  TARGET  DISTANCE 

SCANR,  BEAM  SCAN  R ANGE ( I NCL UD ING  POSITIVE  AND  NEGATIVE  ANGLES) 

CCR  MONITOR  RECEIVER  UPPER  ANGLE  LIMIT 
HTt  ARRAY  TO  STORE  TARGET  HEIGHTS 
ORTt  ARRAY  TO  STORE  TARGET  RANGE 
f MC  t  OPERATING  FREQUENCY  IN  MHZ 
RLL  =  1 1  READ  IN  HEIGHT  AND  DISTANCE 
RLL=2  ,  RANDOM  HEIGHT  AND  DISTANCE 
RLL  =  3t  RE  AO- IN  ANGLES 

IF  THE  FIRST  ANGLE  READING  IS  GREAT  THAN  lOOOt  READS  INCREMENT 
OF  ANGLES.  THE  SECOND  NUMBER  IS  THE  INITIAL  ANGLE  AND  THE  THIRD 
NUMBER  IS  THE  INCREMENT  OF  ANGLES  IN  DEGREES 
RLL  =4  t  RANDOM  ANGLES 

P  TMPP  S  t  POLE  TEST  MULTIPLE  PATH  PHASE 

**.*«*«. ****XF  READ  ANGLE  OR  HEIGHT  AND  RANGEt  ENTER  DATA  CARDT  HERE 
USE  FIC.6  FORMAT 

=  OR  CONVENTIONAL  ARRAYS  t IF  ARRAY  PATTERN  IS  PLOTTED,  ENTER  DATA  FOR 
FINE  STEPS  BEYOND  THE  POSITION  ANGLE  JUST  READ  IN,  USE  15  FORMAT 


............FIFTH  DATA  CARD,  USE  15  FORMAT***************************** 

NCAL ,  NO.  OF  FAILURES  IN  THE  ARRAY 
NCFL ,  NO.  OF  FAILURES  j.N  COMPACT  NETWORK 

****•■»***».. SIXTH  DATA  CARD  USE  F10.6  FORMAT*************************** 

AER ,  ST.  DEVIATION  OF  AMPLITUDE  ERROR  IN  PER  CENT, PHASE  SHIFTER 
PER,  ST. DEVIATION  OF  ANGLE  ERROR  IN  DEGREES,  PHASE  SHIFTER 
?  ME  AN ,  ST.  DEV.  OF  MEAN  INSERTION  PHASE 
P"PEN,  PROBABILITY  PF  OPEN  CIRCUIT  ON  PHASE  SHIFTERS 
PPSSTc,  PROS  OF  A  PHASE  SHIFTER  STUCK  TO  A  CERTAIN  PHASE  STATE 
THE  PROS.  OF  EACH  PHASE  STATE  IS  UNIFORM 
PTNST  FIXED  INSERTION  PHASE  FOR  PHASE  CYCLING 
PEBIT,  ST. DEVIATION  OF  PHASE  ERROR  IN  EACH  PHASE  SHIFTER  BIT 

IF  PEBIT  LT  1.,  THIS  ERROR  IS  IN  PER  CENT  OF  THE  PHASE  SHIFTER 
BIT 

IF  PEBIT  GE  TO  1.  THIS  ERROR  IS  IN  DEGREES 

*••.**..». .SEVENTH  OATA  CARO  USE  F10.6  FORMAT************************** 

NEEDED  ONLY  FOR  COMPACT  ARRAYS 

CSER,  ST.  DEVIATION  OF  AMPLITUDE  IN  PER  CENT,  COUPLE  NETWORK 
CPER,  ST.  DEVIATION  OF  ANGLE  ERROR  IN  DEGREES,  COUPLE  NETWORK 

**.*.*. *•****••****...*.*********•«***•**•**•***************•*•**•**«*. 

IF  HAZETINE'S  STERRING  ALGORITHM  IS  USED  ENTER  THE  SEQUENCE  HERE 
USE  12  FORMAT 


If  Che  program  is  required  Co  compuCe  Che  scacic  array  paccern  or  Co 
ploc  Che  dynamic  paccern  or  Co  generace  Che  phase  shifcer  swicching  histo- 
grara,  Che  program  will  be  cerminaced  afCer  Che  respecCive  Cask,  is  com¬ 
pleted. 

10.  SUMMARY 

In  summary,  chis  simulacion  program  has  Che  following  oucpucs. 

a.  SCaCiscical  discribucion  of  angle  errors  for  carget  at  random 
posicions. 

b.  Statistical  distribution  of  angle  errors  for  fixed  elevation  and 
range  or  at  a  fixed  azimuth  angle  (pole  test). 

c.  Dynamic  scan  pattern  at  a  fixed  angle.  Before  and  after  fil¬ 
tering  for  both  "TO"  and  "FRO"  cases. 

d.  Output  of  a  monitor  which  is  located  at  a  close-in  range  (phase 
errors  are  corrected). 

e.  Array  element  pattern  and  array  pattern  for  both  conventional 
arrays  and  COMPACT  arrays. 

The  condition  in  which  this  simulation  is  performed  is  summarized 
as  follows: 

a.  Type  of  array  (conventional  or  COMPACT). 

b.  Number  of  elements. 

c.  Phase  and  amplitude  error  is  assumed  to  have  a  normal  distribu¬ 
tion  (input  st.  deviation). 

d.  Phase  shifter  truncation  (input  no.  of  phase  shifter  bits). 

e.  Multiple-path  effect  (input  reflection  coefficient  and  phase). 

f.  Band-pass  filter  to  smooth  data  (input  cut-off  frequency  of  the 
filter) . 

g.  Failures  in  phase  shifter  and  in  COMPACT  network  (input  number 
of  failures  and  type  of  failures). 

h.  Variation  of  beam  scan  algorithms.  (Four  algorithms  are 
included:  Bendix,  Hazeltine  No.  1,  Hazeltine  No.  2,  Random.) 

i.  For  pole  test  the  angle  error  is  averaged  over  a  number  of  runs, 
each  with  a  different  constant  insertion  phase.  (Input  No.  of  runs  and 
insertion  phase.) 


j.  Compensation  of  bias  error  of  COMPACT  array  due  to  tilting  of 
element  pattern. 


k.  Each  plot  is  labeled  with  all  the  simulation  conditions. 

l.  Scan  limit. 

11.  REFERENCES: 

1.  "Guidance  Accuracy  Considerations  for  the  Microwave  Landing 
System,"  Robert  J.  Kelly,  NAECON  76,  Dayton,  Ohio.  18  May  1976. 

2.  "Comparison  Study  of  MLS  Airborne  Signal  Processing  Techniques," 
R.  J.  Kelly  and  E.  F.  C.  Laberge,  NAECON  78,  Dayton,  Ohio,  16- 
18  May  1978. 

3.  "MLS  Basic  Wide  Configuration"  Bendix  Communication  Division, 
Proposal  No.  1601,  Jan  4,  1977. 
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Fig.  E-2(a)  —  Multipath  reflection  tor  close  receiver. 
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Fig.  K-3  —  Statistical  distribution  of  ani;lc  errors. 
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Fig.  E-4(a)  —  Pole  test  results. 


Fig.  E-4(b)  —  Pole  test  results  —  bias  error  compensated. 


CONVENT l UNAL  ARKRY-  PHASE  SHIP  TERS-96  3PACING<>v>-'0.60 
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Fig.  E-5  —  Dynamic  pattern  of  a  scanning  array. 


COMPACT  ARRAY:  TILT  ANGLE < DEG.  >rQ.  PHASE  SHIP TERS-24  SPACINGOJ-0. 60 

NO  MULTIPLE  PATH:  OPT.  ERE. (MHZ)- 5060 

COMPACT  NETWORK'  cKAMP.  >-  0.00:  a< PHASE  IN  DEG. )=  0.:  NO.  OP  fAILURES-"  0 
PHASE  SHIT  TER;  a(AMP. 0.00:  u<PHASE  IN  DEG.;-  0.:  FAILURES"  0:  ^bUS 
ANT.  HEIGHT  <PT)-9  RANGE  (FT) -"150 


RNGLE  DEV. (DEGREES) 

Fig.  E-6  —  Monitor  output. 


COMPACT  ARPAY:  TILT  ANGLE (DEG. >-12.  PHASE  SHIP  TERS-24  SPACING<X)r0.60 
NQ  MULTIPLE  PATH-  OPT.  FRE. (MHZ)-506Q 

COMPACT  NETWORK;  cKAMP.  >-  O.QO:  a< PHASE  IN  DEG. >=  0.:  NO.  OF  FAILURES-  0: 
PHASE  SHIFTER:  0(AMP.  )  =  0.00:  oCPHASE  IN  OEG.  )-■  0. :  FAILURES-  0:  4BITS. 


Fig.  E-7  —  Element  pattern  of  a  COMPACT  array. 


NO.  OF  SWITCHED  PHASE  SHIFTERS 

Fig.  E-8  —  Histogram  of  phase  shifter  switching. 
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Ftg.  E-10  —  Flow  diagram. 
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APPENDIX  P 

DIGITAL  FILTER  DESIGN 


A  typical  digital  filter  is  shown  in  Fig.  F-l,  where  T represents  the 
delay  time  in  seconds.  This  filter  can  be  realized  either  by  delay  lines 
or  by  digital  units,  provided  the  input  signal  is  sampled  at  a  rate  satis¬ 
fying  the  Nyquist  criterion.  The  impulse  response  of  such  a  filter  is 

N 

h (t)  =  l  a  6  (t-nT)  (FI) 

n 

n=o 

where  <5(t)  is  an  impulse  function. 

The  transfer  function  in  the  frequency  domain  of  such  a  filter  is 
but  the  Fourier  transform  of  the  impulse  response;  hence, 

N  f 

H(f )  -  Z  a  exp  (-j2im  tt~)  (F2) 

n=o  n  lr 

where  f  »  1/x  represents  the  sampling  rate. 

If  a  2- transform,  such  that 

z  «  exp  ( j  27T  f-)  (F3) 

r 

is  performed  on  the  filter  impulse  response,  one  finds 

N 

H(z)  =  l  a  z"n  (F4) 

_  n 
n=o 

If  feedback  is  provided  as  shown  in  Fig.  F-2,  the  transfer  function 
in  the  z  plane  becomes 

H(z)  =  P(z)/Q(z)  (F5) 


where 

N 

P(z)  -  I  a  z"n  (F5a) 

n 


N 

Q(z)  *  Z  b  a 


(F5b) 


Thus  a  digital  filter  can  be  represented  by  a  rational  function  in 
the  z-domain.  The  design  goal  is  to  find  such  a  rational  function 
which  has  a  prescribed  amplitude  behavior  as  a  function  of  the  frequency 
f/fr.  The  transformation  of  equation  (F3)  establishes  the  relation 
between  the  frequency  f/fr  and  the  variable  z. 

The  transformation  specified  by  equation  (F3)  essentially  maps  the 
entire  frequency  axis  (f)  into  a  unit  circle  in  the  z,  plane.  Further¬ 
more,  the  magnitude  of  the  transfer  function  shown  in  equation  (F5)  has 
a  periodicity  of  fr.  Thus  the  absolute  value  of  equation  (F5)  has  an 
identical  value  at  f  and  at  f  +  nfr  where  n  is  an  integer.  Due  to  this 
periodicity,  the  frequency  response  can  be  viewed  as  limited  to  the 
range , 


0  <  f/f  <1 
—  r  — 

A  digital  filter  in  most  applications  is  required  to  realize  the 
properties  of  a  conventional  bandpass  filter.  An  approach  to  design 
such  a  digital  filter  is  hence  to  transform  a  known  bandpass  filter  to 
a  digital  filter.  Most  conventional  bandpass  filter  characteristics 
are  specified  as  a  function  of  the  imaginary  axis  (a>  axis)  in  the  S 
plane,  which  extend  from  -a>  to  +“  . 

Therefore,  it  is  evident  that  in  order  to  design  a  digital  filter 
having  a  similar  property,  one  is  required  to  find  a  suitable  trans¬ 
formation  which  will  map  the  entire  axis  (co  axis)  in  the  S  plane  into 
a  unit  circle  in  the  z  plane.  Then,  by  use  of  this  transformation,  the 
properties  of  a  bandpass  filter  which  is  synthesized  in  the  S  domain 
can  be  transferred  into  the  z  domain.  Two  pairs  of  such  transformations 
which  can  achieve  such  a  purpose  are  listed  in  the  following: 


(a)  z 


ft  +  S 
ft  -  S 


s 


=  ft 


z  -  1 
z  +  1 


(b)  z 


S  +  ft 
S  -  ft 


s 


=  ft 


z  +  1 
z  -  1 


Where  ft  is  an  arbitrary  constant. 


(F6a) 

(F6b) 


(F7a) 


(F7b) 


We  shall  now  discuss  the  significance  of  these  transformations. 
From  equation  (F6a),  if 


i  =  exp  (j4>) 
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then 


S  *  j  2  tan  ^ 

Thus  the  entire  imaginary  axis  in  the  S  plane  is  mapped  onto  a  unit 
circle  in  the  ?.  plane.  Since  ft  is  an  arbitrary  constant,  one  may  assume 
that 

*c 

ft  =  cot  — 

Under  this  condition,  an  arc  between  +  (j>  on  the  unit  circle  in  the  z 
plane  is  mapped  into  a  region,  S  =  +  j  in  the  S  plane  (see  Fig.  F-3). 

According  to  equation  (F3) ,  frequencies  such  as  iji  >  — p—  correspond  to 

c  —  r 

the  region  in  which  |u)[<_  1.  Thus  a  low-pass  filter  designed  in  the  S 
domain  having  a  cutoff  frequency  of  unity  can  be  transformed  into  the  z 
domain  where  it  will  also  be  a  low-pass  filter  having  a  corresponding 
cutoff  frequency  at 


2irf 


Similarly,  a  high-pass  filter  designed  in  the  S  domain  can  be  transformed 
into  the  z  domain  to  be  a  corresponding  high-pass  filter. 

By  similar  reasoning,  one  finds  that  according  to  equation  (77a) 

_  (J) 

S  *  -j  ft  cot  Y 

If  one  sets 

ft  =  tan  —  » 

one  finds  that  an  arc  between  +  $  such  that  |<f>|  _>  <J>c  will  be  mapped 
into  a  region  S  =  +  j  as  shown  in  Fig.  F-4.  Hence,  in  this  case,  a  low- 
pass  filter  having  a  cut-off  frequency  at  unity  in  the  S  domain  will 
become  a  high-pass  filter  in  the  z domain  with  a  cut-off  frequency 


2irf 


r 


On  the  other  hand  a  high-pass  filter  designed  in  the  S  domain  will 
be  transformed  into  a  low-pass  filter  in  the  z  domain. 

The  above  transformations  furnish  a  convenient  way  to  transer  any 
type  of  bandpass  filter,  either  high-pass  or  low-pass. 
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The  Butterworth  low  pass  filcer  has  a  transfer  function 


T(ju)‘ 


1  +  £  (to/ Cjl)  ) 

C 


2N 


where  N  is  the  order  of  the  rational  function.  The  pole  locations  in 
the  S  plane  are 


S 


k 


* 


where 


k  =  1,2, ...,N 


and  zeros  at  infinity. 


One  may  transfer  these  poles  and  zeros  in  the  S  plane  into  the  z 
plane  by  use  of  equation  (F6a).  From  these  poles  and  zeros  in  the  z 
plane,  a  corresponding  transfer  function  can  be  formed. 
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APPENDIX  G 


Optimum  Beam  Steering  Technique 


1.  Introduction 

This  appendix  describes  a  beam  steering  technique  for  a  linear 
phased  array  antenna  in  which  the  phases  of  the  individual  phase  shifters 
are  determined' by  digital  phase  shifters  with  discrete  settings.  The 
objective  of  the  beam  steering  technique  is  to  provide  maximum  precision 
in  beam  pointing  and  minimum  degradation  in  beam  shape.  The  beam  steer¬ 
ing  technique  described  in  this  appendix  is  felt  to  be  optimum  in  the 
sense  that  it  minimizes  the  energy  lost  from  the  designed  aperture  distri¬ 
bution.  The  technique  uses  a  "recipe"  consisting  of  a  sequence  of 
element  numbers  and  phase  settings  for  smooth  sequential  angular  scan  of 
the  array,  but  the  technique  is  not  limited  to  a  quasi-continuous  scan¬ 
ning  application. 

2.  Theory  of  Operation 

The  antenna  system  under  consideration  is  shown  in  Pig.  G-l.  It 
consists  of  an  array  of  N  elements,  where  N  can  be  even  or  odd.  Each 
element  is  fed  from  the  RF  distribution  network  through  a  digitally  con¬ 
trolled  phase  shifter.  Associated  with  each  element  is  a  quiescent 
excitation,  obtained  when  its  phase  shifter  is  set  at  zero.  The  ampli¬ 
tude  is  determined  by  the  aperture  distribution  associated  with  the 
required  pattern  characteristics.  The  quiescent  phase  may  be  intention¬ 
ally  nonzero  in  the  event  that  "stairstep"  phase  distributions  are  to  be 
avoided.  No  random  errors  are  assumed  in  the  analysis.  It  is  assumed 
that  all  settings  of  all  phase  shifters  are  exact,  and  that  the  overall 
condition  of  the  array  is  always  exactly  calculable. 

A  digital  phase  shifter  is  defined  as  one  which  has  2^  phase  states 
with  a  transmission  phase  difference  between  adjacent  phase  states  of 
2vr/2J  rad.  Thus,  the  transmission  phase  of  a  four-bit  phase  shifter  can 
be  set  to  any  one  of  16  values,  ranging  from  0  to  15tt/ 8  rad  in  tt/8  rad 
steps. 

The  array  can  be  stepped  through  the  phase  states  shown  in  Fig.  G-2 
and  put  into  a  condition  which  is  equivalent  to  the  quiescent  condition. 
Phase  shifters  to  the  right  of  array  center  are  stepped  positively  in 
phase,  and  phase  shifters  to  the  left  are  stepped  negatively  in  phase. 

The  condition  for  which  adjacent  elements  differ  in  phase  by  one  phase 
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step  is  electrically  equivalent  to  the  original  quiescent  condition. 
Thus,  it  is  only  necessary  to  develop  a  beam  steering  technique  for  the 
beam  motion  implied  by  the  phase  states  of  Fig.  G-2.  The  numbers  in  the 
boxes  represent  the  sequence  in  which  the  phase  shifters  are  stepped. 
Symmetrically  opposed  phase  states  are  stepped  alternately. 


If  the  number  of  bits  in  the  phase  shifter  is  J,  and  the  element 
spacing  is  s,  the  beam  motion  shown  in  Fig.  G-2  after  21  phase  steps  is 

A(sin0)  =  *  2^  =  (G-1) 

2J  s2J 

where  9  is  beam  angle  from  Fig.  G-l.  Equation  (G-l)  is  based  on  the 
relationship  for  the  beam  pointing  direction  of  an  array  with  linear 
phase  progression  of  4>  rad  between  adjacent  elements,  which  is 


2tts 

X 


sin0  =  4> . 


Therefore,  the  change  in  sin0  for  a  given  change  in  <)>  is  given  by 


A(sin0)  =  • 

For  the  example  of  Fig.  G-2,  A<J)  =  2tt/2^.  The  maximum  number  of  steps 
that  can  be  used  to  scan  the  beam  through  A(sin0)  is 


«•  (¥)  (¥)-  ^ 


(N  odd) 


Thus,  the  average  step  size  is 

Vsi»e,  -  m&l  .  » 

s2 


(N2-1) 


«  4X 


0JM2 

s2  N 


(N  large) 


(G-2) 


It  is  seen  from  Fig.  G-2  that  the  available  number  of  steps  is  cut 
in  half  by  operating  the  phase  shifters  in  pairs.  It  is  seen  from 
Equation  (G-2)  that  the  step  size  is  doubled  by  removing  one  bit  from  the 
phase  shifter. 


The  design  problem  addressed  in  this  appendix  is  the  choice  of  an 
optimum  technique  for  stepping  the  phase  shifters.  This  technique  is 
based  on  a  single  criterion:  the  total  error  energy  will  always  be  min¬ 
imized.  Error  energy  is  defined  as  the  sum  of  the  squares  of  the  com¬ 
ponents  obtained  by  taking  the  difference  between  actual  element  excita¬ 
tion  and  a  weighted  LI1S  fitted  linear  phase  function.  In  addition,  the 
beam  direction  is  defined  by  the  linear  phase  function. 
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To  derive  the  steering  technique,  we  use  y  as  the  phase  coordinate 
and  x  as  the  element  location  coordinate.  The  location  of  the  ith  ele¬ 
ment  and  its  associated  phase  are  x^  and  y^,  respectively.  The  amplitude 
excitation  is  a^,  and  the  power  excitation  is  w^  =  . 

Although  we  allow  for  independent  control  of  all  phase  shifters, 
we  claim  that  the  preferred  phase  distribution  is  always  symmetric,  and 
we  calculate  only  symmetric  distributions.  Symmetric  pairs  can  still 
be  shifted  in  sequence. 

We  first  derive  a  linear  phase  function,  y  =  bx,  such  that  the 
weighted  sum  of  the  squared  error  components  is  minimum. 

The  phase  error  for  the  ith  component  is 


Ay, 


y.  -  bx. 
i  1 


The  energy  in  the  error  component  of  the  ith  element  is 


E.  =  (a. Ay.)2 
i  ii 

•  =  w^ (Ayi ) 2 

=  Vyk2“2bxiyi+b2xi2) 

The  total  energy  is  to  be  minimized  with  respect  to  b: 
E  *  Zwi(y.2-2bxiyi+b2xi2) 

||  =*  Zw.(-2xiyi+2bxi2)  S*C  0 
2 

bZw. x.  =  Ew. x.y. 
i  i  ill 


(G-3) 


b 


SwiVi 


Ew.x. 

i  i 


(G-4) 


Substitution  of  Equation  (G-4)  into  Equation  (G-3)  yields 


If  we  now  change  the  phase  of  the  jth  element  by  6y.»  the  change 


in  E  is  <5Ej  ,  and  is  given  by 


J 


5E,  =  w.6y. 
J  J  J 


Ew.x.y . 
111 

yj'Xjl - 2 

EVl 


w.x. 

+  5y_.  <*  1 - 3  -j 

Ew.x 
i  i 


(G-5) 


Equation  (G-5)  allows  us  to  apply  the  minimum  energy  criterion  by 
always  selecting  that  j  which  gives  the  minimum  value  of  <5E..  In  addi¬ 
tion,  it  is  found  that  beam  motion  is  given  by  3 

<Sy.  w.x. 

6b  =  3  33  (G-6) 

J  Ew.x. 

l  l 

Note  that  this  beam  steering  technique  is  deterministic.  There  is 
always  a  preferred  phase  shifter  to  step  (sometimes  more  than  one  with 
the  same  minimum  6Ej ) .  Thus,  once  an  initial  set  of  y^  and  an  aperture 
distribution  w-^  are  specified,  the  beam  stepping  sequence  is  determined. 

Table  G-I  gives  samples  of  beam  stepping  sequences  for  three  cases 
for  a  13-element  array.  For  the  first  case,  the  amplitude  distribution 
is  uniform  (w^  =  1)  and  the  initial  phases  are  zero  (y^  =  0).  For  the 
second  case  the  amplitude  distribution  has  been  tapered  to  a  cosine- 
squared  on  a  pedestal  in  power.  The  values  of  w^  are  determined  by 

w*  *  .3  +  .7  cos-2  (~TT~)  •  For  the  third  case,  the  taper  is  retained, 
and  a  quasi-random  set  of  initial  phases  has  been  introduced.  The 
values  of  yt  for  i  =  1  through  6  are  .42,  -  .25,  .08,  .25,  -  .42,  and 
-  .08.  The  value  of  6y^  is  taken  to  be  unity,  so  that  one  complete  scan 
cycle  is  contained  between  b  =  0  and  b  =  1.  The  beam  stepping  sequence 
is  thereby  normalized  with  respect  to  the  number  of  bits  in  the  phase 
shifters.  That  is,  the  same  sequence  is  used  regardless  of  the  number 
of  bits.  The  beam  simply  doesn't  move  as  far  if  there  are  more  bits. 


Note  that,  for  a  given  aperture  distribution,  the  size  of  the  beam 
step  is  determined  by  the  element  that  is  shifted  as  predicted  by  Equa¬ 
tion  (G-6).  Note  also  that,  for  the  cases  for  which  all  y^  are  initial¬ 
ly  zero,  the  stepping  sequence  is  essentially  symmetrical  about  the  11th 
step.  On  the  10th,  11th,  and  12th  steps,  elements  1,  3,  and  5  must  be 
shifted,  and  no  symmetry  is  possible  for  that  portion  of  the  sequence. 


It  is  readily  seen  that  the  stepping  sequence  is  strongly  dependent 
on  the  initial  values  of  y^  and  less  strongly  dependent  on  the  aperture 
distribution.  It  is  also  pointed  out  that,  whereas  the  ratio  of  maximum 
to  minimum  beam  step  sizes  is  (N-l)/2  (for  N  odd)  or  N-l  (for  N  even), 
for  a  uniform  aperture  distribution,  it  is  always  less  for  a  tapered 
aperture  distribution. 
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Table  G-I.  Examples  of  Beam  Steering  Technique 
Applied  to  a  13-Element  Array 


Step; 

-i 

j 

Uniform  Amplitude 

7i  =  0 
b  AE 

1  J 

Taper 
=  0 
b 

AE 

y± 

i 

Taper 
+  0;  b0  = 
b 

-.0073 

AE 

i 

i 

6 

.0659 

.6044 

'  6 

.0467 

.2232  ! 

i 

5 

.0413 

-.0038 

2 

5 

.1209 

.0659 

5 

.0953 

.1125  ‘ 

6 

.0880 

.0199 

3 

4 

.1648 

-.1429 

■  4 

.1481 

.0138 

2 

.1306 

.0534 

4 

3 

.1978 

-.0879 

;  3 

.2002 

-.0311 

4 

.1834 

.1284 

5 

2 

.2198 

.1648 

i  2 

.2428 

.0968 

3 

.2355 

-.0668 

6 

6 

.2857 

-.0330 

6 

.2895 

-.0600 

5 

.2841 

-.1699 

7 

5 

.3407 

-.1319 

'  .  5 

.3381 

-.0535 

6 

.3308 

-.2633 

8 

4 

.3846 

.0989 

4 

.3909 

,0441 

6 

.3775 

.1831 

9 

6 

.4505 

-.0110 

•  6 

.4376 

.0090 

5 

.4261 

.0552 

10 

1 

.4615 

.0879 

1 

.4617 

.0967 

4 

.4789 

-.0632 

11 

5 

.5165 

.1099 

'  3 

.5137 

.0510 

3 

.5310 

.0902 

12 

3 

.5495 

-.1978 

5 

.5624 

-.1478 

6 

.5777 

.0583 

13 

6 

.6154 

.0110 

:  6 

.6091 

-.0090 

5 

.6263 

.0544 

14 

4 

.6593 

-.0989 

4 

.6619 

-.0441 

2 

.6689 

-.0768 

15 

5 

.7143 

.1319 

5 

.7105 

.0535 

4 

.7217 

-.0329 

16  ' 

6 

.7802 

.0330 

6 

.7572 

.0600 

6 

.7684 

-.0311 

17  ; 

2 

.8022 

-.1648 

2 

.7998 

-.0968 

5 

.8170 

.0904 

18 

3 

.8352 

.0879 

3 

.8519 

.0311 

3 

.8691 

.0703 

19 

4 

.8791 

.1429 

4 

.9047 

-.0138 

6 

.9158 

.0405 

20 

5 

.9341 

-.0659 

5 

.9533 

-.1125  : 

4 

.9686 

-.0197 

21 

6 

1.0000 

-.6044 

6 

1.000 

-.2232  • 

1 

.9927 

-.1164 
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3.  Hardware  Application 

In  order  to  establish  a  quasi-continuous  scanning  pattern,  a  func¬ 
tional  relationship  must  be  defined  between  b  and  time,  and  the  appro¬ 
priate  phase  shifters  will  be  shifted  at  the  resulting  times.  In  this 
arrangement  both  time  and  element  values  must  be  stored  in  the  beam 
steering  prescription.  In  practice  it  might  be  preferable  to  execute 
phase  shifts  at  uniform  time  intervals,  corresponding  to  a  clock  rate. 

For  example,  the  13-element  array  depicted  in  Table  G-I  actually  takes 
42  steps  to  scan  from  b  -  0  to  b  =  1,  counting  both  positive  and  negative 
element  positions.  If  we  specify  that  we  want  some  number  of  evenly 
spaced  steps,  it  is  possible  to  combine  groups  of  phase  shifts  in  order 
to  closely  approximate  the  desired  beam  motion.  There  are  a  number  of 
procedures  by  which  steps  could  be  combined,  but  the  most  straightforward 
appears  to  be  to  take  the  tabulated  steps  in  order  and  group  them  so  that 
the  actual  beam  position  during  a  particular  dwell  period  is  as  close  to 
the  desired  position  as  possible.  For  example,  if  we  want  ten  beam 
positions  between  b  =  0  and  b  =  1  uniformly  spaced  in  time.  Figure.  G-3 
shows  that  the  desired  setting  values  of  b  should  be  .05,  .15,  .25,... 
.95.  For  the  uniform  amplitude  and  taper  with  y^  £  0  cases  of  Table  G-I, 
the  resulting  process  is  shown  in  Table  G-II.  The  error  in  b  is  shown 
for  each  beam  position,  and  the  RMS  error  in  b  is  .0172  for  the  uniform 
case  and  .0162  for  the  tapered  case.  Based  on  21  total  available  steps, 
which  will  give  a  maximum  deviation  from  any  desired  setting  of  about 
1/42  ■  .0238,  and  assuming  a  uniform  distribution  over  +  .0238,  we 
would  expect  an  RMS  error  of  .0238/v^  =  .0137,  which  is  in  reasonable 
agreement  with  the  results  obtained. 

So  long  as  the  number  of  steps  used  is  less  than  the  total  number 
available,  we  would  expect  the  RMS  error  to  be  independent  of  the  number 
used.  That  is,  the  RMS  error  in  sin8  is  approximately 

a  ■  X/  (s2J+1M/3)  *  A/  (s2J-1S2v/3) 

where  we  have  inserted  the  step  size  from  Equations  (G-l)  and  (G-2). 

Since  the  beamwidth  B  in  sin9  is  approximately  A/Ns,  the  precision  of 
beam  steering  relative  to  a  beamwidth  is 

ff/B  -  1/2J"1N^. 

Two  likely  hardware  realizations  of  these  techniques  both  require 
the  memorization  of  the  stepping  sequence  in  some  form — it  is  not  reason¬ 
able  to  expect  to  compute  the  sequences  during  the  scan.  One  procedure, 
which  is  more  applicable  to  a  small  number  of  phase  shifters  and/or  a 
limited  angular  scan,  is  to  commit  the  entire  scanning  sequence  to  memory 
and  use  it  as  a  recipe  from  a  ROM.  Another  procedure  is  to  commit  a 
single  portion  of  the  repeating  sequence  to  memory  and  to  calculate  suc¬ 
cessive  values  of  b,  using  the  memorized  sequence  to  determine  the  phase 
shifter  steps  for  each  clock  interval.  A  third  possibility  is  becoming 


more  feasible  as  the  cost  of  semiconductor  memory  continues  to  decline. 
The  entire  stepping  sequence  for  given  phase  shifter  can  be  located  at 
the  phase  shifter  so  that  the  central  beam  steering  unit  would  consist 
of  little  more  than  an  initializer  and  a  clock. 

Table  G-II.  Scanning  Sequences  for 
Uniform  Step  Size 


Desired 

b 

Uniform  Amplitude 
Realized  Units 

b  Shifted 

Error 

Desired 

b 

Taper,  b  = 
Realized 
b 

-  0.0073 
Units 
Shifted 

Error 

.05 

.0659 

6,6 

+.0159 

.043 

.0413 

1,5 

.15 

.1648 

4,5 

+.0148 

.143 

.1306 

2,6 

-.0124 

.25 

.2198 

2,3 

.243 

.2355 

3,4 

.35 

.3407 

5,6 

-.0093 

.343 

.3308 

5,6 

-.0122 

.45 

.4505 

4,6 

+.0005 

.443 

.4261 

5,6 

-.0169 

.55 

.5495 

1,3,5 

•  .543 

.5310 

3,4 

-.0120 

.65 

.6593 

4,6 

+.0093 

.643 

.6263 

5,6 

-.0167 

.75 

.7802 

5,6 

+.0302 

.743 

.7217 

2,4 

-.0213 

.85 

.8352 

2,3 

-.0148 

.843 

.8691 

3,5,6 

+.0261 

.95 

.9341 

4,5 

-.0159 

.943 

.9686 

4,6 

+.0256 

For  the  recipe  approach,  three  examples  are  given,  ranging  from 
limited  coverage  and  small  number  of  phase  shifters  to  full  coverage  and 
large  array.  Table  G-III  lists  the  parameters  and  the  resulting  required 
memory  capacity. 

Table  G-III.  Memory  Required  for  Recipe  Method  of  Beam 
Steering  Control;  4-bit  Phase  Shifters 

No.  Phase  Angle  Memory  Size 

Shifters  Coverage  (deg)  (k  bits) 


30 

0-15 

2.7 

100 

0-15 

35 

100 

-40  -  +40 

180 

76 


For  the  second  technique,  which  is  part  ROM,  part  processing,  the 
memory  requirement  for  the  ROM  is  about  8  k  bits  for  the  100-element 
array  of  4-bit  phase  shifters  covering  +  40  degrees. 

For  the  case  in  which  all  memory  is  located  at  the  phase  shifters, 
the  memory  requirement  at  each  phase  shifter  would  be  about  3  k  bits  for 
the  100  element  array  of  4-bit  phase  shifters  covering  +  40  degrees. 
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APPENDIX  H 

Design  of  a  360-Degree  Azimuth  Antenna 
For  the  Microwave  Landing  System 


lo  Introduction 


The  purpose  of  this  appendix  is  to  outline  the  theoretical  and  design 
considerations  involved  in  the  achievement  of  an  electronically  scanned 
360-degree-coverage  azimuth  antenna. 

The  appendix  is  presented  in  three  sections.  The  first  section  after 
this  one  considers  the  basic  array  geometry  and  concludes  that  a  circular 
array  of  radiating  elements  is  preferable  to  a  number  of  linear  arrays. 
The  next  section  addresses  several  alternative  scanning  techniques  for 
circular  arrays.  The  final  section  describes  the  particular  scanning 
circular  array  that  is  being  developed  at  NRL.  Appendices  I  and  J  pre¬ 
sent  the  theory  of  a  compact  3D  lens  configuration  which  will  be  used  in 
the  electronic  scanning  feed  system. 

It  is  assumed  that  the  antenna  will  satisfy  the  following  system 
requirements : 


Azimuth  beamwidth  4  degrees 

Azimuth  sidelobe  level  25  dE 

Elevation  pattern  Shaped 

Theoretical  gain  31  dB 

Antenna  Losses  4  dB 

Net  maximum  gain  27  dB 


2„  Array  Geometry 

Two  basic  array  geometries  are  available  for  achieving  360-degree 
azimuth  coverage  with  an  electronically  scanned  antenna — a  circle  or  a 
polygon.  The  circular  array  is  a  single  array  with  a  single  scanning 
feed  system.  The  polygon  is  a  set  of  independent  linear  arrays. 

From  the  standpoint  of  system  performance,  the  two  geometries  have 
•  ;<f  following  features: 

«  ~he  beamwidth  changes  with  azimuth  angle  for  the  linear  arrays, 
•'tmwidth  of  the  circular  array  degrades  slightly  with  elevation 

r  tills  application,  for  example,  the  phase  error  at  15  degrees 
.  -i  -28  degrees. 
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h  The  circular  array  produces  a  "planar"  beam— that  is,  beam 

varies  slightly  with  elevation  angle. 

The  conical  beams  produced  by  the  linear  arrays  will  produce  an 
.rrof  in  tSe  azimuch  tadlcatad  by  MIS  recovers.  Over  an  azlautb  region 
which  increases  with  Increasing  elevation  angle,  the  “'borne  tece  ^ 
will  receive  two  mainbeam  passages  rather  than  one. 
receiver  ignores  the  second  beam,  an  error  is  generated,  the  maxi 
value  of  which  increases  rapidly  with  elevation  angle.  The maximum 
error  for  two  elevation  angles,  together  with  the  azimuth  angle 
which  the  error  is  obtained,  are  listed  below.  It  is  assumed  that  four 
linear  arrays  are  used  and  that  they  point  in  the  azimuth  directions,  n90 


Elevation 

Azimuth 

Error 

20° 

n90°-45°+3.8° 

.42° 

30° 

n90°-45°+9. 7° 

2.24° 

The  use  of  linear  arrays 

implies  that  a 

number  of  independent  elec- 

tronically  scanned  arrays  will  be  used.  The  total  number  of  array  ele- 

ments,  and  consequently  phase 
beamvddth. 

shifters,  is 

listed  below  for  a  four-degree 

Number  of 

Elements  per 

Total  Number 

Arrays 

Array 

of  Elements 

3 

60 

180 

4 

39 

156 

5 

32 

160 

6 

28 

168 

It  is  seen  that  the  case  of  four  arrays  minimizes  the  number  of 
phase  shifters.  However,  the  corresponding  calculation  for  a  circular 
array  results  in  an  array  of  about  108  elements.  It  will  be  shown  that 
the  most  complex  beam  scanning  system  which  could  be  used  with  the  cir¬ 
cular  array  would  require  about  94  phase  shifters.  The  selected  circu¬ 
lar  array  scanning  configuration  requires  37  phase  shifters  and  37  3-way 
switches. 

Therefore,  the  use  of  linear  arrays  involves  both  errors  and  increased 
hardware  complexity,  and  the  remainder  of  this  appendix  is  devoted  to  the 
design  of  a  scanning  circular  array. 

It  is  assumed  that  the  circular  array  design  will  be  based  on  the 
following  criteria: 
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a. 

ture. 


Aperture/beamwidth.  9.  ,  (deg)  =  64A/A,  where  A  is  effective  aper- 
This  relationship  is  determined  by  the  25  dB  sidelobe  requirement. 


b,  Aperture/diameter.  A  =  D  cos  9  ,  where  D  is  diameter  of  the  cir¬ 
cular  array  and  9  is  half  the  effective  coverage  angle  of  the  circular 
array  elements. 

c.  Element  spacing.  S  =  — ^  q-  ■  .  This  relation  comes  directly 

from  array  theorv,  where  9  is  defined  above. 

-  »  e 


The  above  criteria,  together  with  the  assumed  system  requirements  and 
with  9g  =  60°,  result  in  the  following  antenna  parameters: 


93dB 


=  4° 


Beamwidth 


A 

D 

S 


N 


16A  Effective  aperture 

18. 5A  Array  diameter 
.536A  Element  spacing 
108  Number  of  elements 


Dimensions  1.1  m.  diameter,  ~  1  m. 


height 


3.  Alternative  Circular-Array  Scanning  Techniques 

Two  basic  electronic-scan  techniques  are  considered  here.  One  is 
directly  analogous  to  the  phase-scanned  linear  array  and  will  be  referred 
to  as  phase  scan.  We  define  this  as  a  configuration  in  which  the  indi¬ 
vidual  element  patterns  are  identical  in  amplitude  but  are  orthongonal 
because  of  their  differing  phase  centers.  The  other,  which  we  will 
refer  to  as  amplitude  scan,  is  directly  analogous  to  InterScan.  In  this 
configuration,  the  patterns  are  orthogonal  based  on  their  amplitude 
characteristics  alone. 


The  radiating  elements  of  a  circular  array  do  not  satisfy  either  of 
the  above  criteria.  The  patterns  are  orthogonal,  but  not  based  on  either 
phase  or  amplitude  alone.  Thus,  there  exists  an  essentially  continuous 
range  of  orthogonal  patterns  available  from  the  circular  array,  with 
phase  orthogonal  at  one  extreme  and  amplitude  orthogonal  at  the  other. 

We  will  examine  techniques  for  transforming  from  the  element  patterns  to 
any  other  set  of  orthogonal  patterns. 

Returning  to  the  two  types  of  electronic  scan  techniques,  it  is 
noted  that  the  formation  of  pattern  functions  by  these  two  methods  repre¬ 
sents  the  two  classical  synthesis  procedures.  In  phase  scan,  the  pat¬ 
tern  is  synthesized  from  element  patterns  which  can  be  described  as  com¬ 
plex  Fourier  harmonics,  and  this  is  a  Fourier  synthesis.  In  amplitude 
scan,  the  pattern  is  synthesized  in  a  Woodward  fashion  from  element 
patterns  which  are  essentially  sin  x/x  functions.  In  the  case  of  the 
Fourier  synthesis,  the  number  of  harmonics  used  depends  on  the  size  of 
the  aperture.  In  the  Woodward  synthesis,  the  number  of  sin  x/x  patterns 
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used  depends  on  the  precision  required  of  the  synthesis.  An  advantage 
of  amplitude  scan  is  that  it  requires  fewer  degrees  of  freedom  for  the 
synthesis  of  a  directive  pattern  with  control  over  a  few  sidelobes. 

As  is  well  known  to  workers  in  the  field,  amplitude  scan  consists  of 
a  coarse  scan  mechanism  and  a  fine  scan  modulator.  The  coarse  scan  is 
accomplished  by  a  switching  network  which  selects  the  patterns  that  will 
be  used  in  the  synthesis.  The  fine  scan  modulator  selects  the  specific 
amplitude  weights  that  will  be  applied  to  the  patterns.  Now,  in  order 
to  use  phase  control  rather  than  amplitude  control  and  thereby  afford 
straightforward  comparison  among  the  alternative  techniques,  we  choose 
to  Fourier  transform  the  amplitude  weighting  function  supplied  by  the 
fine-scan  modulator.  This  weighting  function  is  sampled  according  to 
the  pattern  locations,  and  the  sampling  points  translate  with  the  beam 
scanning  motion.  The  Fourier  transform  changes  translation  of  one 
function  into  phase  modulation  of  another  function. 

One  might  ask  why  we  have  bothered  to  use  an  amplitude  scan  config¬ 
uration  if  we  are  going  to  transform  it  into  a  phase  scan  arrangement. 

The  reason  for  using  amplitude  scan  is  the  reduction  in  complexity. 

Fewer  phase  control  devices  will  be  required  for  the  amplitude  scanned 
antenna. 

Scanning  configurations  are  now  described  for  the  circular  array. 
Figure  H-l  lists  the  four  patterns  which  will  be  used  as  examples.  The 
patterns  range  from  omnidirectional  to  the  minimum  beamwidth  patterns 
available  from  the  array,  and  the  scanning  techniques  range  from  pure 
phase  scan  to  pure  amplitude  scan.  The  second  technique  uses  the  pattern 
of  the  circular  array  radiating  element  directly,  which  has  a  beamwidth 
of  about  90  degrees  and  useful  coverage  of  about  120  degrees.  All  of 
the  other  patterns  are  obtained  by  suitably  transforming  the  array  ele¬ 
ment  patterns. 

For  the  pure  phase  scan  case,  the  use  of  omnidirectional  patterns 
represents  the  transformation  of  the  circular  array  to  an  equivalent 
linear  array,  and  Figure  H-2  indicates  how  these  patterns  would  be  gen¬ 
erated.  A  mode  forming  network  such  as  a  Butler  matrix  or  a  Rotman  lens 
is  used. 

For  the  amplitude  scan  case,  a  microwave  optical  collimator  or  an 
appropriate  network  must  be  used  to  form  the  multiple  directive  patterns. 
Possible  techniques  are  the  Myers  geodesic  lens,  the  Luneberg  lens,  or 
transformation  from  the  equivalent  linear  array  to  multiple  patterns.. 

The  last  example  involves  two  transformations,  one  to  the  equivalent  V 
linear  array  and  another  to  the  multiple  beams.  Each  transformation 
involves  a  Fourier  transforming  network  such  as  a  Rotman  lens  or  a  Butler 
matrix.  Figure  H-3  outlines  the  resulting  configuration,  showing  the 
multiple-beam-forming  device,  the  coarse  switching  network,  and  the  fine 
scan  modulator,  with  the  number  of  interconnections.  Using  still  another 


but  small  Fourier  transform  network,  the  fine  scan  modulator  is  realized 
with  four  phase  shifters.  Therefore,  it  is  seen  that  phase  scan  uses  a 
maximum  number  of  phase  shifters  and  no  switches  while  amplitude  scan 
uses  a  minimum  number  of  phase  shifters  and  a  maximum  number  of  switches. 

Examples  of  two  intermediate  techniques  are  given.  One  uses  the 
element  pattern  directly  without  translation  and  requires  moderate  numbers 
of  phase  shifters  and  switches.  This  is  referred  to  as  a  cyclic  scan 
technique.  The  other  uses  slight  transformation  of  the  element  patterns, 
also  requires  moderate  number  of  phase  shifters  and  switches,  and  is 
called  a  hybrid  scan  technique. 

The  hybrid  scan  technique  uses  ten  phase  shifters  and  transformed 
patterns  with  about  36  degrees  beamwidth.  These  patterns  can  be  obtained 
by  means  of  a  partial  focus  lens,  analogous  to  the  Luneberg  lens,  which 
does  not  achieve  perfect  focus.  The  detailed  characteristics  of  this 
lens  are  not  known.  However,  the  same  performance  can  be  obtained  with 
an  equivalent  linear  array  and  a  set  of  small  Fourier  transform  networks. 
Although  the  mix  of  switches  and  phase  shifters  for  this  configuration 
might  be  desirable,  the  lack  of  an  available,  practical  device  to  achieve 
the  partial  focus  precludes  its  consideration  at  this  time. 

The  pure  phase  scan  configuration,  although  it  requires  no  switches, 
involves  the  use  of  a  large  network  or  lens  to  provide  the  equivalent 
linear  array  outputs.  In  addition,  it  requires  about  100  phase  shifters. 
This  degree  of  complexity  allows  complete  pattern  synthesis — in  theory, 
every  sidelobe  of  the  360-degree  pattern  can  be  specified.  Because  of 
the  excessive  lens  size  and  number  of  control  devices,  this  alternative 
is  rejected. 

Figure  H-4  illustrates  the  details  of  the  various  scanning  tech¬ 
niques,  assuming  particular  realizations  for  the  lenses  and  transforming 
networks.  The  remaining  choice  lies  between  amplitude  or  InterScan  and 
cyclic  scan.  The  choice  is  somewhat  subjective.  If  we  prefer  reduced 
sensitivity  to  component  failure,  we  choose  cyclic  scan  with  36  separate 
switches  and  36  separate  phase  shifters.  This  alternative  is  selected 
and  is  described  in  the  following  section. 

4.  Cyclic  Scan  Circular  Array 

Referring  again  to  Figure  H-4,  it  is  noted  that  cyclic  scan  is  func¬ 
tionally  similar  to  the  other  techniques  except  for  the  lack  of  a  trans¬ 
forming  lens  between  the  array  and  the  coarse  scan  switches.  Furthermore, 
the  fine  scan  network  is  the  major  part  of  the  configuration,  consisting 
of  a  lens,  a  power  divider,  and  36  phase  shifters. 

As  was  pointed  out  above,  a  smaller  number  of  control  devices  pro¬ 
duces  a  lesser  degree  of  control  over  the  synthesized  pattern.  The  use 
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of  36  phase  shifters  will  provide  control  over  about  36  sidelobes,  or 
about  120  degrees  of  the  360  degree  antenna  pattern.  This  is  completely 
adequate  to  provide  the  desired  beamwidth  and  sidelobe  level. 


Figure  H-4  indicates  a  Rotinan  lens  with  dimensions  larger  than  the 
array  for  the  cyclic  scan  system.  The  Rotman  lens  is  a  Fourier  trans¬ 
forming  device  which  is  planar  in  construction,  generally  using  a 
parallel-plate  transmission  medium.  An  alternative  three-dimensional  lens 
has  been  devised,  which  can  provide  the  same  transfer  characteristics  as 
the  Rotman  lens  in  a  more  compact  manner.  This  lens  concept  is  described 
in  Appendices  I  and  J0  It  is  shown  that  the  three-dimensional  lens  must 
have  37  ports,  and  its  dimensions  are  about  .6  m,  which  is  significantly 
smaller  than  the  circular  array.  Thus,  the  entire  feed  system  can  be 
packaged  within  Che  circular  array  aperture. 

The  resulting  antenna  design  has  the  following  parameters: 


Circular  array  diameter 

Number  of  elements 

Element  spacing 

Switches 

Phase  shifters 

Lens 


1.09  m. 

Ill 

3.09  cm. 

37,  1P3T 
37 

Cylindrical,  ~.6  m.  long  by 
~,6  m.  diameter.  37  ports 
on  each  side 


i 


108  CABLES 

93-  108  CABLES 

4*  1:23-27 
4  CABLES 

4  PHASE  SHIFTERS 


Fig.  H-3  —  InterScan  technique. 
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Fig.  H-4  —  Comparison  of  scanning  techniques. 


APPENDIX  I 

THREE-DIMENSIONAL  BOOTLACE  LENSES 


Introduction.  Previously  published  work,  on  bootlace  lenses  has 
predominantly  treated  two-dimensional  configurations  and  linear 
antenna  arrays.  A  three-dimensional  bootlace  lens  offers  the 
principal  advantage  of  achieving  multiple  beams  from  a  planar 
array  with  frequency-independent  collimation  and  beam  positions. 

When  used  with  circularly  polarized  lens  and  feed  radiators,  small 
phase  adjustments  are  easily  nade  by  rotating  the  radiators. 

Because  of  the  inherent  polarization  diversity  of  the  transmission 
medium,  the  use  of  dual-polarized  lens  and  feed  radiators  enables 
the  realization  of  two  lenses  in  one  structure,  capable  of  inde¬ 
pendently  feeding  two  arrays . 

Two  different  design  approaches  to  the  three-dimensional  bootlace 
lens  are  presented  and  analyzed  here.  Equations  and  plots  of  per¬ 
formance  and  design  parameters  are  presented.  A  technique  is 
described  in  Appendix  J  for  appropriately  configuring  a  3D  bootlace 
lens  so  that  its  outputs  can  be  connected  to  a  linear  array  to 
produce  multiple  beams.  The  result  of  this  procedure  is  to  reduce 
greatly  the  overall  size  of  lens  required  to  obtain  multiple 
beams  from  a  linear  array. 

Air-Filled  3D  Lens.  It  is  well  known  that  an  air-filled  3D  boot¬ 
lace  lens  exhibits  astigmatism  is  its  dominant  aberration  for  off- 
axis  feed  positions.  Although  it  is  possible  to  achieve  a  quadru- 
focal  design  with  one  focus  on-axis  and  the  other  three  symmetri¬ 
cally  spaced  off-axis,  it  is  found  that  minimum  aberration  is 
obtained  with  all  foci  collocated  on  the  lens  axis. 

A  "universal"  3D  air-filled  bootlace  lens  design  has  been  derived 
by  requiring  that,  for  any  radial  distance  off-axis,  with  both 
feed  point  and  lens  surface  points  located  at  this  radius,  the 
astigmatism  should  be  minimum.  The  lens  and  feed  surfaces  are 
also  required  to  be  symmetric.  The  cross-section  of  the  result¬ 
ing  lens  is  shown  in  Figure  1-1.  Points  on  the  lens  surface  are 
defined  by  (x,z),  and  points  on  the  feed  surface  are  defined  by 
(x,-z).  The  axial  points  of  these  surfaces  are  +z  .  The  lens 
surface  is  generated  by  the  equation 

(x2+z2)^  +  z  +  (2x2+4z2)^  =  4 [x2+(zq+z)2] 2  -4zq  (1-1) 

The  length  of  the  transmission  line  between  a  point  on  the  lens 
surface  and  the  aperture  is  given  by 

L (x)  -  2zq  -  [x2+  (zo+z)2]^  (1-2) 
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Fig.  1-2  —  Characteristics  of  3D  air-filled  bootlace  lenses. 


If  Che  lens  diameter  is  taken  to  be  a  particular  value  of  x,  then 
the  aperture  size  for  hemispheric  coverage  is 

A(xx)  =  2[(x2+z2)2  -z.J  (1-3) 

In  the  plots  of  Figure  1-2,  all  parameters  are  scaled  for  unity 
aperture  diameter  and  hemispheric  coverage.  The  astigmatic  wave- 
front  error  is  cosinusoidal  in  <5  and  exhibits  a  maximum  deviation 
of  +E  relative  to  the  aperture  size.  In  the  limit  of  large  D,  E 
varies  inversely  as  D2  .  For  the  lens  configurations  of  most 
likely  practical  interest,  E  ranges  from  .1  to  .01. 

3D  Fisheye  Lens .  The  lens  configuration  of  Figure  1-3  has  inhomo¬ 
geneous  dielectric  loading  in  the  region  between  the  feed  and 
lens  surfaces.  This  region  is  obtained  by  cutting  out  a  portion 
of  a  Maxwell  fisheye  lens,  which  has  a  refractive  index  given  by 
n  =  2/(l+r2).[l]  The  lens  and  feed  surfaces  are  spherical  with 
radii  of  /I,  defined  by 

r2  +  2rcos  <J>  -1  =  0 

All  transmission  lines  are  of  equal  length,  and  the  lens  has  one 
feed  point  of  perfect  focus  on  the  axis.  It  can  be  shown  that 
the  lens  is  anastigmatic,  and  the  dominant  off-axis  aberration 
is  coma..  The  lens  ana  feed  surfaces  are  geodesics — that  is,  the 
surfaces  contain  ray  paths.  Thus,  the  surfaces  appear  to  be 
planar  to  the  radiators  located  on  them.  Furthermore,  all  rays 
through  the  center  of  one  surface  intersect  the  opposite  surface 
normally.  In  this  regard,  the  lens  is  optically  equivalent  to 
the  2D  spherical  geodesic  bootlace  lens  proposed  by  Wild  [2]. 

Analysis  of  the  optical  performance  of  this  lens  requires  the 
calculation  of  the  optical  path  length  between  any  two  points  in 
the  Maxwell  fisheye  medium.  In  general,  it  is  very  difficul.  to 
calculate  optical  path  lengths  in  inhomogeneous  media.  However, 
because  the  fisheye  is  an  ideal  optical  system  and  because  all 
ray  paths  are  circles,  it  is  always  possible  to  find  an  arc  on 
the  unit  circle  that  is  optically  equal  to  a  given  ray  path.  In 
this  way  the  lens  can  be  analyzed  by  straightforward  geometric 
techniques.  The  optical  path  length  between  any  two  points  is 
given  by 

Y  =  tt  -  tan  1  [cota  Vl+a2-]  -  tan  ^cotCB-ci)  Vl+a2]  (1-4) 
tana  =  (Rz-R1cos  9)/R1sin9 

a2  =  4 (R2+R2  -  2RiRzcos  0)/sinz  0 

Rj  = 

R2  =  (l-r2)/r2 

where  ^  and  r2  are  the  radial  distances  of  the  two  points  and  9 
is  the  angle  between  them. 
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If  the  lens  and  feed  points  are  defined  by  <p  and  6,  where  <p  is 
shown  in  Figure  1-3  and  6  is  the  angle  between  the  radii  as  seen 
vertically  along  the  axis  (separation  is  minimum  for  5=0  and  max¬ 
imum  for  6=180°),  the  optical  path  length  is  given  by 


tan(y/2) 


l-sin2({icos26/2  ‘ 

2 

J_  l-sin2(J)sin  6/2  J 
The  desired  optical  path  length  is  given  by 


Y  =  tt/2-(L/2)cos  6 
o 


L  *  Tr-4tan  1  (cos 


(1-5) 


(1-6) 


Equations  (1-5)  and  (1-6)  are  used  to  determine  wavefront  error, 
£=Y-Y0. 


As  in  the  case  of  the  3D  air-filled  bootlace  lens,  the  fisheye 
bootlace  lens  is  also  universal  in  that  a  smaller  portion  of  the 
lens  can  be  used  to  obtain  improved  optical  performance.  If  the 
edge  of  the  feed  and  lens  surface  is  defined  by  <j>  and  all  parame¬ 
ters  are  scaled  so  that  the  aperture  diameter  is  unity  and  cover¬ 
age  is  hemispheric,  the  plots  of  Figure  1-4  are  obtained.  Mote 
that  for  <t  <  90°,  the  maximum  value  of  r  is  less  than  one,  so 
that  the  minimum  value  of  the  refractive  index  n  is  greater  than 
one.  The  index  and  lens  dimensions  have  also  been  scaled  in  Fig¬ 
ure  1-4  so  that  the  minimum  index  is  one  and  the  maximum  index  at 
lens  center  is  as  shown.  The  comatic  wavefront  error  is  essen¬ 
tially  cubic  in  cos6.  In  the  limit  of  small  <J)  or  large  D,  E 
varies  inversely  as  .  For  the  lens  configurations  of  most 
likely  practical  interest,  E  ranges  from  .07  to  .0007. 


The  performance  of  the  fisheye  lens,  although  inferior  to  that  of 
the  conventional  2D  bootlace  lens,  is  far  superior  to  that  of  the 
3D  air-filled  lens. [3] 
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APPENDIX  J 

Using  3D  Lenses  to  Obtain  Multiple  Beams  from  Linear  Arrays 

In  Ref.  [1]  Shelton  discussed  the  conditions  under  which  multiple 
beams  can  be  obtained  from  planar  arrays  with  lossless  hybrid  networks. 

It  was  shown  that  such  multibeam  systems  are  realizable  only  if  the 
array  satisfies  a  tiling  condition — that  is,  the  infinite  regular  lattice 
of  radiating  element  locations  can  be  completely  filled  by  translations 
of  the  planar  array.  It  was  also  shown  that  certain  planar  arrays  can 
be  fed  by  linear-array  multibeam  networks.  In  particular,  all  square 
arrays  and  all  symmetric  hexagonal  arrays  fall  in  this  special  category. 
Figure  J-l  shows  a  19-element  hexagonal  array  with  adjacent  arrays  in  a 
tiling  configuration.  The  elements  are  numbered  so  that  arithmetic 
modulo-19  progressions  are  obtained  along  any  linear  cut  through  the  array. 
The  modulo-19  sequence  is  constrained  to  the  integers  between  -9  and  +9, 
inclusive.  Figure  J-2  shows  the  beam  positions  of  the  multiple  beams 
from  the  array  of  Figure  J-l  and  some  of  the  grating  lobe  locations  for 
these  beams. 

It  was  shown  in  [1]  that  the  integer  steps  from  element  to  element 
can  be  determined  by  requiring  that  the  tiled  arrays  all  have  the  same 
Integer  patterns  on  them.  If  a  is  the  integer  step  along  the  ax  lattice 
vector  and  b  along  the  a2  lattice  vector,  the  equations  can  be  shown  to 
be 
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Fig.  J'l  —  19-element  hexagonal  array  showing  adjacent  arrays  in  tiling  arrangement. 
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Fig.  J-2  —  Beam  positions  for  array  of  Figure  J-l. 
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5a  -  3b  -  19m 
3a  +  2b  *  19n 


(J-l) 


for  m  =  0,  n  »  1,  equations  (J-l)  yield  a  *  3,  b  »  5. 

These  conditions  can  also  be  applied  to  planar  and  linear  arrays 
fed  by  bootlace  lenses  because  bootlace  lenses  can  be  theoretically 
equivalent  to  a  hybrid  multibeam  network  [2].  The  logical  extension  of 
the  above  conclusions  is  that,  for  the  special  square  and  hexagonal 
arrays,  if  a  three-dimensional  lens  can  generate  multiple  beams,  then 
the  same  lens  can  be  used  to  feed  a  linear  array  with  the  same  number  of 
elements. 

We  use  the  19-element  array  of  Figure  J-l  as  an  example,  and  we 

assume  that  a  3D  bootlace  lens  is  used.  The  feed  and  lens  surfaces  will 

consist  of  19-element  hexagonal  arrays.  The  parameters  which  must  be 

determined  at  this  point  are  the  orientation  of  the  lens  and  feed  arrays, 

the  correct  scaling  of  the  lens  dimensions,  and  the  assignment  of  lens 

and  feed  elements.  The  required  orientation  is  first  derived. 

»  n 

The  grating  lobe  locations  0,  0  ,  and  0  in  Figure  J-2  are  determined 
by  the  geometry  of  the  planar  array  lattice.  The  lattice  vectors  for 
this  array  are  given  by 


sx 


a,  -  s(x/2  +  yi(3/2). 


(J-2) 


where  s  is  spacing  between  nearest  radiators  and  x  and  y  are  unit  vectors 
in  the  x  and  y  directions.  All  element  locations  in  the  lattice  are 
given  by  linear  combinations  of  at  and  a2  : 

re  “  nat  +  ma2 

*  s[(n+  l/2)x  +  myv^3/2]. 


97 


where  ?  is  the  vector  from  the  origin  to  the  element  location  and  m,  n 
assume  all  positive  and  negative  integer  values.  The  grating  lobe  loca¬ 
tions  are  defined  in  the  u  ,  u  plane  by 

x  y 

R  -  R  +  2iTk(u  +  u Jfi)  +  Jlu  (4tt/  /3)  (J-3) 

go  x  Y  y 

where  u  *  (27rs/X)sin6  and  u  •  (2iTs/X)sin0  ,  u  and  u  are  unit  vectors, 
x  x  y  y  x  y 

k,2.  also  assume  all  integer  values,  and  Rq  is  the  location  of  any  one  of 

the  grating  lobe  maxima.  Note  that  the  minimum  nonzero  magnitude  of 

R  -I  is 
g  o 

It  is  clear  that  the  on-axis  beam  is  obtained  with  the  on-axis  feed 


element  and  that  the  off-axis  beams  are  obtained  by  elements  of  the  feed 
array  that  are  correspondingly  displaced  off-axis.  Thus,  if  the  beam 
is  located  at  R  in  the  u  ,  u  plane,  the  feed  element  will  be  located 

s  x  y 

-y 

at  a  point  in  the  projected  feed  array  plane,  x,  y,  given  by  -  c  R^» 
where  c  is  a  constant  of  proportionality. 

Given  the  on-axis  location  of  one  of  the  beams  and  the  lattice  geome¬ 
try  defined  by  equations  02),  the  grating  lobe  maxima  must  be  located 
as  given  by  equation  (J-3)  for  Rq  -  0  and  as  shown  in  Figure  J-2.  The 
locations  of  the  remaining  18  beams  generated  by  the  lens  can  only  be  as 
shown  in  Figure  J-2  (or  as  a  mirror  image  in  the  uy  axis  obtained  by  mul¬ 
tiplying  all  values  of  ux  by  -1). 

The  basis  vectors  for  the  multiple  beam  locations  are  given  by 

r  »  (2iT(2n+l)/N)tt  -  (2tt/ /3N)u  (J_4) 

gi  x  y 

R  -  (2tt (n+l)/N)u  +  (2TT(3n+l)//3N)u  , 

g2  x  y 

for  the  symmetrical  hexagonal  array  with  n  rings  of  elements  and  number 
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of  elements  N  *  3n2+3n+l.  Note  that  |R  I  =  |R  I  *  47t/^5n  •  R  ■ -R  hifa. 

gi  '  g2  go 

As  pointed  out  previously,  the  feed  element  locations  will  be  opposite  the 
multiple  beam  locations.  Whereas  the  radiating  array  and  lens  array  loca¬ 
tions  given  by  equations  (J-2)  lie  on  a  triangular  lattice  with  basis  vec¬ 
tors  oriented  at  integral  multiples  of  60  deg  relative  to  the  x-axis,  the 
feed  array  basis  vectors  given  by  equations  (J-4)  are  rotated  with  respect 
to  the  x-axis  by  an  angle  y  given  by 
coty  *  -  v;3(2n+l), 

for  the  symmetrical  hexagonal  array.  For  the  19-element  array,  y  * 

6.5868  deg. 

It  is  interesting  to  consider  that  for  any  orientation  of  feed  and 
lens  arrays  other  than  the  one  derived  here,  the  antenna  system  will  not 
be  optimum  in  the  sense  of  producing  lossless,  orthogonal,  maximum  gain 
multiple  beams  from  a  planar  array. 

A  similar  analysis  can  be  carried  out  for  arrays  with  n2  elements, 
for  which  the  lens  and  feed  arrays  are  square  and  located  on  square 
element  lattices.  Here  again  it  is  found  that  the  feed  array  must  be 
rotated  with  respect  to  the  lens  array.  In  this  case  the  rotation  angle 
y  is  given  by 

tany  ■  m/n, 

where  m  is  an  integer  constrained  by  0<m<n/2,  and  m  cannot  be  a  factor 
of  n  other  than  one. 

We  now  have  two  steps  remaining  to  complete  the  analysis  of  the  use 
of  3D  lenses  to  form  multiple  beams  from  linear  arrays.  The  first  is 
the  proper  scaling  of  the  lens  to  achieve  the  necessary  phase  relation¬ 
ships.  The  second  is  to  relate  the  lens  elements  to  the  elements  of  the 


99 


linear  array  and  to  relate  the  feed  elements  to  the  beam  positions 
formed  by  the  3D  lens/linear  array  combination. 

Design  formulas  for  3D  lenses  are  based  on  the  assumptions  that 
the  aperture  size  is  unity,  the  maximum  path  delay  across  the  aperture 
for  an  edge  feed  position  is  also  unity,  and  the  feed  and  lens  surfaces 
are  geometrically  similar.  Referring  to  the  modulo-19  number  of  the 
elements  of  the  array  of  Figure  J-l,  we  note  that  these  numbers  should 
be  multiplied  by  ,\/19  to  give  path  delay,  and  that  the  modulo-19  num¬ 
bering  should  first  be  discarded  at  the  edge  of  the  array,  as  shown  in 
Figure  j-3.  The  feed  position  indicated  by  dashed  circle  c  produces  a 
path  delay  of  20A/19  between  lens  elements  a  and  b,  and  the  feed  posi¬ 
tion  indicated  by  dashed  circle  d  produces  a  path  delay  of  40A/19.  The 
geometry  of  the  elements  at  a,  b,  and  d  before  scaling  is  shown  in  Fig¬ 
ure  J-4.  The  gap  between  the  edges  of  the  feed  and  lens  surfaces  is  G, 
the  diameter  of  these  surfaces  is  D,  and  the  rotation  angle  between  lens 
and  feed  arrays  is  y.  The  path  delay  between  paths  ad  and  bd  is  given 
by 


Ail  -  Vg2+D2cos2y/2  -  ^G2+D2sin2Y/2 
We  also  have  Vg2+D2  -  G  -  1  or  D2  -  2G+1. 

It  is  seen  that  Ail  *  1  for  Y  "  0  and  decreases  with  increasing  y.  There¬ 
fore,  we  scale  the  lens  dimensions  by  the  factor  40A/19Ail.  For  lens 
designs  of  maximum  compactness  A£  becomes  significant:  For  G  *  0  and 
Y  *  6.5868  deg,  Ail  *  .9409.  However,  for  a  lens  with  G  -  1,  D  ■  fi. 

Ail  -  .9926. 
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Fig.  J-3  —  Lens  array  showingrelative  locations  of 
elements  a,b  and  feed  points  c,d. 
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For  the  19-element  array  shown  in  Figure  J-l,  the  numbers  show  the 
connections  to  the  elements  of  the  linear  array.  Figure  J-2  shows  the 
positions  of  the  multiple  beams  produced  by  the  linear  array.  If  one 
off-axis  beam  position  can  be  determined,  as  in  the  case  of  beam  number 
one  in  Figure  J-3,  then  all  beam  positions  are  defined  by  extending  the 
linear  progression.  This  results  first  in  numbering  the  row  of  beams 
consisting  of  0  through  -5  in  Figure  J-4,  then  the  rows  -4  to  -8  and  -5 
to  -7,  based  on  the  tiling  characteristics  of  the  grating  lobes  of  the 
multiple  beams.  Continuing,  we  can  number  the  rows  -8  to  7,  8  to  3,  and 
finally  3  to  0'  and  3  to  0”,  to  complete  the  numbers  of  all  beam  posi¬ 
tions.  Note  that  for  both  arrays  the  elements  around  the  center  element 
can  be  defined  by  three  integers,  2,  3,  5  for  the  lens  array  and  1,  7,  8 
for  the  feed  array.  These  sets  of  integers  are  sufficient  to  define  all 
elements  for  the  arrays.  Other  sets  of  integers  can  be  obtained  by 
forming  Integer  multiples,  modulo  19,  of  the  initial  sets.  The  result 
Is  that  lens  and  feed  arrays  can  be  formed  with  any  combination  of  three 
sets  of  integers,  2,  3,  5,  and  1,  7,  8,  and  4,  6,  9. 

This  appendix  has  outlined  a  procedure  for  feeding  a  linear  array 
with  the  outputs  of  a  3D  bootlace  lens.  The  procedure  involves  the  fol¬ 
lowing  steps: 

1.  Select  either  a  square  or  regular  hexagonal  array. 

2.  Determine  the  rotation  angle  between  feed  and  lens  arrays. 

3.  Compute  the  appropriate  scale  factor  for  the  lens. 

4.  Determine  the  connections  from  the  lens  arrays  to  the  linear 
arrays  and  the  beam  positions  assigned  to  elements  of  the  feed  array. 
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